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ABSTRACT 


Since  a  nonlinear  feedback  control  system  may  possess  more  than  one 
type  of  forced  oscillations,  it  is  highly  desirable  to  investigate  the 
type  of  forced  oscillations  which  can  occur  when  the  nonlinear  restoring 
force  function  is  of  a  specific  form. 

This  paper  is  to  present  a  harmonic  linearization  method  for  finding 
the  existence  of  forced  oscillations  and  response  curve  characteristics 
of  a  nonlinear  feedback  control  system  by  means  of  finding  the  restoring 
force  function  of  the  nonlinearity  and  using  the  harmonic  balance  and  an 
iteration  method  for  investigating  the  conditions  for  one  type  of  forced 
oscillations  exhibited. 

The  existence  conditions  for  fundamental  frequency,  sub-harmonic  of 
2nd  order  and  3rd  order  forced  oscillations  of  a  second  order  feedback 
control  system  are  investigated;  also  the  fundamental  frequency  forced 
oscillation  for  a  higher  order  system  and  the  jump  resonance  frequencies 
of  response  curve  are  investigated,  and  a  general  expression  for  the  equi¬ 
valent  gain  of  the  nonlinearity  has  been  developed. 

The  author  wishes  to  express  his  appreciation  for  the  assistance  and 
encouragement  given  by  Dr.  George  Julius  Thaler  of  the  U.  S.  Naval  Post¬ 
graduate  School  in  this  investigation. 
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CHAPTER  I 


INTRODUCTION 

One  of  the  major  branches  of  servo  mechanism  analysis  for  nonlinear 
control  systems  is  the  study  and  methods  of  prevention  of  continuous  oscil¬ 
lation.  For  a  linear  control  system,  the  output  may  oscillate  with  con¬ 
tinuous  or  increasing  amplitude  with  no  input  signal  applied,  in  which  case, 
the  system  is  said  to  be  unstable,  or  they  may  have  oscillations  which  die 
away.  In  both  cases  this  transient  response  is  the  characteristic  of  the 
system  itself  and  the  conditions  for  divergent,  continuous  or  decaying 
oscillation  do  not  depend  on  the  form  or  magnitude  of  the  input  signal. 

When  a  nonlinear  element  is  present,  this  independence  of  the  input 
signal  r.o  longer  holds.  The  characteristic  performance  of  the  nonlinear 
control  system  will  depend  on  both,  i.e.,  forcing  function  and  the 
characteristic  of  the  nonlinear  element. 

A  nonlinear  control  system,  under  suitable  conditions  may  exhibit 
steady  oscillations  in  which  the  main  component  has  a  frequency  which  is 
dependent  on  the  frequency  of  the  forcing  function.  This  type  of  oscilla¬ 
tion  is  called  "Forced  Oscillation".  A  forced  oscillation  for  which  the 
frequency  is  a  fraction  of  the  forcing  function  frequency  is  called  "Sub¬ 
harmonic  Forced  Oscillation".  A  forced  oscillation  for  which  the  frequency 
is  multiple  of  the  forcing  frequency  is  called  "Super  Harmonic  Forced  Oscil 
lation". 

Since  a  nonlinear  control  system  may  possess  more  than  one  type  of 
forced  oscillations,  it  is  highly  desirable  to  investigate  the  type  of 
forced  oscillations  which  can  occur  when  the  nonlinear  restoring  force 
function  is  of  a  specific  form. 

Some  methods  of  investigation  for  nonlinear  forced  oscillation  have 
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been  discussed  by  J.  C.  West1  in  his  book  "Analytical  Techniques  for 

Nonlinear  Control  System"  (1960).  A  graphic  method  for  investigating 

forced  oscillation  of  "ON  OFF"  control  system  has  been  developed  by  Hamel 

2 

(French)  and  Tsypkin  (USSR)  .  Some  work  for  subharmonic  oscillations  of 

a  typical  restoring  force  function  has  been  done  by  C.  A.  Ludeke  and 

3  „  4 

William  Pong  in  their  paper  1959*  and  also  by  Ogata  in  his  PhD  thesis  , 

June  1956. 

There  are  some  other  authors  for  investigation  of  the  vibration  of  a 

5  6  7 

mechanical  system,  such  as  N.  Minorsky  ,  J.  J.  Stoker  ,  and  Y.  H„  Ku  in 
their  books  of  "Introduction  to  Nonlinear  Mechanics" ;"Nonlinear  Vibrations" 
and  "Analysis  and  Control  of  Nonlinear  System"  respectively. 

The  purpose  of  this  paper  is  to  present  a  harmonic  linearization 

g 

method  for  finding  the  existence  of  forced  oscillations  and  response 

curve  characteristics  of  a  nonlinear  feedback  control  system,  by  means  of 

9  10 

finding  the  restoring  force  function  of  the  nonlinearity  9  ,  and  using  an 

iteration  method  for  investigating  the  condition  for  one  type  of  forced 
oscillations  exhibited.  A  response  curve  and  "Jump  Resonance"  can  be 
also  investigated. 

The  emphasis  of  this  paper  is  placed  on  showing  the  general  approach 
of  the  method,  therefore*  second  order  nonlinear  feedback  control  systems 
are  used  as  examples.  However*  a  higher  order  system  will  be  discussed  in 
this  paper  also. 

By  the  harmonic  linearization  and  iteration  method*  a  general  equa¬ 
tion  for  the  equivalent  gain  of  an  odd  function  nonlinear  element  will  be 
developed* *.  It  is  in  terms  of  the  amplitude  and  the  frequency  of  input. 

There  are  six  chapters  in  the  main  body  of  this  paper.  The  first 
chapter  and  the  last  chapter  are  a  general  introduction  and  conclusion 
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respectively.  Chapter  two  is  a  general  description  of  oscillation  of 
feedback  control  systems  and  an  investigation  of  the  restoring  force  func¬ 
tion  of  the  nonlinear  element.  Chapter  three  investigates  the  existence  of 
fundamental  forced  oscillations  and  the  jump  phenomena  in  the  nonlinear 
feedback  control  system.  Chapter  four  is  a  sub  harmonic  oscillation  in¬ 
vestigation.  Chapter  five  investigates  the  existence  of  forced  oscillation 
of  a  higher  order  feedback  control  system. 


a 
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CHAPTER  II 


GENERAL  DISCUSSION  OF  OSCILLATION  IN  FEEDBACK  CONTROL  SYSTEMS 
2-1  General  Description: 

For  either  a  linear  or  nonlinear  feedback  control  system,  most  re- 
quirements  for  control  system  design  specify  a  minimum  time  response,  high 
accuracy,  and  high  stability.  It  is  clearly  in  order  to  have  a  minimum 
time  response  and  high  accuracy  for  a  feedback  control  system,  but  first 
of  all,  the  system  should  be  stable. 

In  general,  two  types  of  oscillations  may  be  exhibited  in  a  feedback 
control  system,  one  is  free  oscillation,  i.e.,  an  autonomous  system.  This 
is  a  system  which  is  left  by  itself  with  its  outside  source  of  energy  sudden 
ly  cut  off  or  abruptly  changed  in  amplitude.  In  such  a  system,  free  oscilla 
tions  occur  as  the  system  tends  to  equilibrium  after  reaching  a  state  of  un¬ 
equilibrium.  The  oscillation  frequency  of  this  type  depends  only  on  the 
characteristic  of  system  itself.  Ordinarily,  a  free  oscillation  is  a  damp¬ 
ed  oscillation;  as  the  time  increases  the  amplitude  decreases  and  damps  out 
in  steady  state.  For  an  ideal  case,  If  it  is  a  system  without  damping,  i.e. 
oC  »  0,  the  oscillation  will  be  continuous. 

The  other  type  of  oscillation  is  a  forced  oscillation,  i.e.  a  nonauto- 
nomous  system.  This  is  a  system  that  is  acted  upon  by  an  outside  source  of 
energy,  a  force  function.  The  force  function  may  be  a  constant;  a  periodic 
function  or  any  other  function  of  time. 

The  analysis  and  methods  for  preventing  oscillations  of  a  feedback 
control  system  are  concerned  directly  or  indirectly  with  the  characteristic 
of  system,  i.e.,  the  differential  equations  of  the  system.  If  the  system 
is  linear  a  frequency  response  of  the  system  can  be  directly  solved  from 

f 

the  differential  equation,  furthermore,  the  principle  of  superposition  can 
be  applied  too.  On  the  other  hand,  if  the  system  is  nonlinear,  either  with 
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a  nonlinear  damping  or  nonlinear  restoring  force  function*  then  the  basic 
tools  for  the  analysis  of  a  linear  system  are  no  longer  valid, 

2-2  Oscillation  in  Linear  Systems: 

There  are  two  types  of  methods  for  analysis  and  design  of  linear  feed¬ 
back  control  systems;  one  is  graphic  method,  the  other  is  an  analytical 
method.  The  best  known  of  graphic  methods  consist  of  Bode  diagram;  Polar 
plot;  Nichols  diagram;  and  Root  locus.  If  we  know  the  transfer  function  of 
the  system,  any  of  the  above  methods  can  be  used  for  solving  the  problem  of 
stability.  The  analytical  method  is  a  mathematical  analysis  for  solving  the 
differential  equation,  and  plot  the  curves  of  the  solved  equation. 

First  consider  a  simple  linear  feedback  control  system,  as  shown  in 
Fig.  2-1.  The  forward  element  has  a  transfer  function  G(s),  and  a  feedback 


O.iS) 


CrCSJ 


r 

i 


1  h  > ) 

L _ 


Fig.  2-1  Block  Diagram  of  a  Linear  System 
element  transfer  function  H(s).  The  feedback  signal  is  H(s)0c  (s)  and  it 
is  assumed  that  this  is  subtracted  from  the  input  signal  (s)  to  form 

the  control  signal  e(s);  thus: 

e(s)  S  0r( s)  -  H(s)  £>c  (s)  (2*1) 

The  relation  between  the  control  signal  and  reference  signal  in  the 
Laplace  transformation  form  is: 

4  <*>  -  —  £$<;■(,) - ^<S)  (2‘2) 
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The  solution  of  equation  (2-2)  for  a  particular  input  can  be  obtained 
by  Laplace  and  Heaviside  methods.  This  involves  the  determination  of  the 
roots  of  equation: 

1  +  H(s)G(s)  =  0  (2-3) 


In  order  to  split  the  transfer  function  into  the  sum  of  partial  fractions 
of  standard  form  for  which  time  solutions  are  known.  If  equation  (2-3)  is 
an  n  th  order  polynominal  in  s,  there  will  be  n  roots  say  to  such 

that  (s-  S  )  is  a  factor  of  the  polynominal.  Special  cases  arise  in  which 
some  of  roots  are  identical,  and  if  m  roots  of  occurs,  (s  -  S t)m  is  a 
factor  of  the  polynominal. 

In  general,  the  roots  will  be  a  complex  number 8  ass 


S  »  +  j  cu 

r  r  J  r 


(2-4) 


in  which  U*  and  are  both  real  numbers,  and  since  the  coefficients  of 

r  r 

the  polynominal  of  equation  (2-3)  are  all  necessarily  real,  then  the  com¬ 
plex  roots  should  be  in  conjugate  pairs,  Thus: 


(s  *  +  j  «Jr)  (  S  -  *c  r 


j  CO  ) 

J  r/ 


will  be  a  factor  of  the  polynominal. 

In  determining  the  transient  response  to  a  forcing  function  input, 
these  various  factors  become  the  denominator  of  the  partial  fraction  re¬ 
presentation  of  the  transfer  function  equation  (2-2).  The  term  of  K^/ (s  -S^) 
gives  a  time  solution  containing  expoential  term: 

e(^  r  +  J  co  r)t 

The  complete  solution  is  the  sum  of  all  such  terms. 

The  complex  conjugate  pairs  give  rise  to  an  oscillatory  term: 

e  *  Sin(<y  t  +  ) 

If  the  system  has  convergent  response,  all  of  these  terms  in  the  time  solu¬ 
tion  must  converge  and  hence  all  roots  S,  to  S  must  have  negative  real 

1  n 

parts.  If  any  of  these  roots  has  a  zero  real  part,  then  a  continuous 
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oscillation  is  produced*  If  any  of  the  complex  roots  has  a  positive  real 
part,  then  divergent  oscillation  will  be  the  result* 

To  determine  the  stability  of  a  system  by  using  the  graphic  method,  it 
is  not  necessary  to  find  the  roots  of  the  characteristic  equation*  If  the 
equation  is  known  algebrically,  the  Routh  method  can  be  applied  and  also 
by  the  Nyquist  and  Bode  diagram  method  a  graphical  solution  can  be  obtained* 
Nyquist  criteron  graphical  analysis  is  the  most  useful  method  to  solve 
the  stability  problem  of  control  systems*  Rewriting  equation  (2-3)  in  the 
form  of: 

G(s)H(s)  — 1  (2-5) 

The  roots  of  equation  of  (2-3)  are  particular  values  of  s  which  satisfy 
equation  (2-5).  Thus  if  the  equation  of  G(s)H(s)  can  be  mapped  for  continu¬ 
ously  variable  values  of  s,  then  those  values  for  which  G(s)H(s)  ■  -1  can 
be  determined. 

Assume  a  particular  roots  of  +  j  can  be  represented  by  a 

particular  value  of  G(s)H(s),  i.e. 

G(^r  +  J  f)H(.i  r  +  j^r)  -  R(4.  &•■■)  +  j  I<*'r  .  i  )  (2>6) 

This  equation  represents  a  vector  in  I  vs  R  plane,  the  magnitude  of  vector: 

I  2  2 

Ms  Jr  +  I  (2-7) 

and  the  direction  by  an  angle: 

*  Tan"1l/R  (2-8) 

For  a  constant  ,  the  value  of  M  will  be  a  locus  in  I~R  plane  as 
Co  varies.  This  process  can  be  repeated  for  different  values  of <6  r  to 
obtain  a  family  of  curves  as  shown  in  Fig.  2-2* 

This  steady  state  frequency  response  locus  of  the  open  loop  system 
can  be  obtained  experimentally  from  the  system  by  measuring  gain  and  phase 
for  a  sinusodial  input  of  varying  frequency*  It  can  be  shown  that  is  the 
steady  state  frequency  response  locus  passes  to  the  right  side  of  (-1*  0) 
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Fig*  2-2  Locus  of  an  Open  Loop  Transfer 
Function  G(s)H(s) 

point,  this  system  will  be  stable,  on  the  other  hand,  if  it  passes 
through  or  to  the  left  side  of  (-1,  0)  point,  the  system  will  be  unstable* 
This  criterion  is  due  to  Nyquist.  A  typical  curve  is  shown  in  Fig*  2-3. 

In  the  extension  to  nonlinear  control  systems,  or  more  complex  multi¬ 
loop  systems,  the  single  loop  criterion  is  sufficient  for  the  majority  of 
needs  where  G(s)  can  represent  the  over  all  transfer  function  of  several 
loops  provided  that  G(s)  is  itself  stable. 


Fig.  2-3  Nyquist  Criterion  for  Three 
Characteristic  States 
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2-3  Forced  Oscillations  in  a  Linear  System: 

Assume  a  second  order  linear  system  block  diagram  is  shown  in  Fig* 

2-1 ;  in  which: 

g<s>  -  . s  (K-—  (2-9) 

and: 

H(s)  -  1  (2-10) 

From  equation  (2-2) t  the  differential  equation  of  system  becomes: 

0  +  0  +  K9  *  K0  (2-11) 

c  c  c  r 

Let  the  forcing  function: 

0  (t)  *  F  Cos  to  t  (2-12) 

Equation  (2-11)  becomes: 

9  +  X  0  +  K0  *  C  Cos  U)  t  (2-13) 

c  c  c 

in  which  C  is  a  constant  of  value  of  FK« 

The  solution  of  equation  (2-13)  consists  of  the  sum  of  the  solution 
of  the  homogeneous  equation*  (i«e.  the  free  oscillation  of  the  system)  and 
the  solution  of  the  non-homogeneous  equation*  Let  the  solution  of  homo¬ 
geneous  equation: 

ec(t)  *  e“2  t(C1  CosiO  ./  1-/  t  +  C2  Sin  it'n/  l“/s  t)  (2-14) 

in  which: 

(On  •  K 

cJ_  m  2;'V  ’L 

f  *  Cos 

The  complete  solution  of  equation  (2-13)  will  be: 

e  (t)  *  e  -- ‘(C.  Cos  a-  :jlf  t  +  C  Sinw;  /lTt)  +— 7===^====?  (2- 15) 

Equation  (2-15)  is  obtained  by  a  superposition  of  the  free  oscillation 
and  the  forced  oscillation  which  varies  from  the  action  of  external  force: 

The  frequency  of  forced  oscillation  is  the  same  as  that  of  the  external 
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(2-16) 


force,  the  magnitude  of  the  forced  oscillation  is  given  by: 

H-  . C - 

(  l  *V7  J  1  ^  J  ^  7  ^ 

The  value  of  phase  shift  8  related  to  the  external  force  is  given 

by: 

Cose  -  - - — --------  -  (2-17) 

,  ?  -  H/  )"  f  x  f  '  h-  * 

Sine  -  — - - —  (2-18) 

;  ■'  if  'riS  hi  ■ 

In  the  case  of  positive  damping,  i.e.  U  y  0,  it  is  clear  from  equation 
(2-15)  that  after  a  sufficiently  long  period  of  time,  the  free  oscillation 
is  damped  out  and  only  the  forced  oscillation  would  be  observed* 

In  the  case  of  no  damping,  i.e.  u  ®  0,  the  phase  shift  0  is  seen 
from  equation  (2-17)  and  (2-18)  to  be  zero  for  OJ  <  ^  and  7t  for  cOs-Xv  In 
other  words,  the  forced  oscillation  is  in  phase  with  the  external  force  if 
the  forced  oscillation  frequency  is  less  than  the  free  oscillation  fre¬ 
quency,  and  is  180  degrees  out  of  phase  with  the  external  force  when  u)  is 
greater  than  . 

In  the  case  of  ^  *  0,  equation  (2-15)  becomes: 

0  (t)  *  (C.  CosfV  t  +  C0  SintO  t) 
c  l  n  2  n 

+  C  Cos  ( cl)  t  +  8)  (2-19) 

\XX  ~  H)*\ 

There  will  be  a  superposition  oscillation  of  two  frequencies,  if  the  value 
of  co  X-  «/S  one  of  which  is  the  natural  frequency,  and  the  other  is  the 
frequency  of  external  force.  In  the  case  of  u)  * tQ *  ,  the  free  and  the  forced 
oscillations  have  the  same  frequency,  and  the  solution  of  equation  (2-13) 
will  be  found: 

C 

0  (t)  *  C-  Cosc^t  +  C0  S  in  uj  t  -  ^;t  Sintdt  (2-20) 

cl  2  2& 

In  this  case,  the  component  due  to  the  external  force  is  no  longer  periodic, 
it  is  oscillatory  with  an  amplitude  that  increases  linearly  with  time.  At 
this  condition,  it  is  a  resonance  phenomenon.  It  is  very  important  in  de¬ 
sign  to  avoid  this  phenomenon  in  a  system  with  a  periodic  force  function. 


10 


In  the  case  of  a  damped  system,  i.e.,  u  >  0,  from  equation  (2-16)  for 
the  steady  state,  the  amplitude  of  forced  oscillation  is  always  finite®  It 
is  possible  to  give,  by  the  use  of  dimensionless  variable,  a  more  general 
significance  to  the  expression  for  the  amplitude  of  forced  oscillation,  let 
the  amplitude: 

C 

|H|  -  M(  --  )  (2-21) 

Therefore  equation  (2-16)  becomes: 

M  -  7~==r - — - r  (2-22) 

v. '  1  rK p*  ( Z~)L 

Where  M  is  defined  to  be  the  magnification  factor®  The  extreme  values  for 

2  2  2 

M  are  attained  for  tO  *  0  and  (  co  f  tC%  )  *1-2  jJ  .  If  1  -  2  f  «<.  0, 

2 

there  is  a  maximum  for  cO  *  0;  if  1  -  2  j°  >  0,  and  0,  there  is  a  maxi¬ 
mum  for  us/cVh  =  /V  '  and  a  minimum  for  co  *  0.  For  small  values  of  damp¬ 

ing  coefficient,  the  frequency  which  produced  the  maximum  amplitude  is  very 
nearly  the  natural  frequency  of  the  system.  Fig®  2-4  shows  the  forced  oscil¬ 
lation  response  curves  of  a  linear  system  as  a  function  of  t-u/  u)n  with 
various  values  of 

2-4  Oscillations  in  Nonlinear  Feedback  Control  System: 

As  mentioned  in  section  2-1  of  this  chapter,  when  a  nonlinearity  is 
present  in  a  feedback  control  system,  the  properties  of  proportionality 
and  superposition  are  not  valid.  In  the  linear  system,  these  properties 
involve  the  existence  of  a  transfer  function  and  of  characteristic  fre¬ 
quencies  proper  to  the  system  (frequencies  at  which  the  system  tends  to 
oscillate  with  an  amplitude  depending  on  initial  conditions);  in  the  case 
of  nonlinear  systems,  the  amplitude,  like  frequency,  can  depend  at  one  and 
the  same  time  on  both  the  initial  conditions  and  the  system  itself.  Some¬ 
times,  as  in  the  case  of  limit  cycle  both  frequency  and  amplitude  of  input 
are  characteristic  of  the  system  and  independent  of  of  initial  conditions. 
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Curve  For  a  Lino 
C'-ViJJ/xiion.s 


In  many  cases  the  analysis  and  design  of  such  systems  using  linear 
theory  may  produce  an  excellent  result.  This  general  procedure  is  a 
type  of  linearization  procedure,  i.e«,  the  actual  nonlinear  system  is  re¬ 
placed  by  a  linear  system,  which  approximates  it,  and  the  analysis  and  de¬ 
sign  techniques  are  applied  to  the  linear  equivalent.  The  most  useful  method 
for  nonlinear  linearization  is  that  of  the  "Describing  Function"  method.  It 
is  actually  a  first  harmonic  approximation  method. 

The  method  for  linearization  used  by  this  paper  is  a  multi-harmonic 
linearization  and  iteration  method;  first  assume  a  harmonic  solution  for 
the  system  differential  equation  and  insert  it  in  the  differential  equation, 
and  compare  the  coefficients  of  the  same  order  harmonic  terms0  A  detailed 
procedure  will  be  discussed  in  the  next  chapter. 

As  previously  explained,  a  nonlinear  feedback  control  system  may  posses 
some  type  of  oscillation.  Which  type  will  exist  depends  on  the  restoring 
force  function  of  nonlinear  element  and  the  amplitude  of  forcing  function. 
Hence,  before  investigating  the  existence  of  forced  oscillations,  the  re¬ 
storing  forcing  function  of  the  nonlinearity  should  be  investigated. 

2-5  Characteristic  of  Nonlinear  Elements: 

Nonlinear  elements  which  may  be  seen  in  feedback  control  systems  can  be 
represented  by  a  combination  of  four  fundamental  concepts  of  nonlinearity, 
which  are:  (1)  Saturation;  (2)  Variable  gain;  (3)  Dead  zone;  (4)  Hysteresis. 
These  basic  notions  of  nonlinearity  can,  in  combination  with  one  another, 
result  in  almost  all  types  of  nonlinearity  in  practice,  as  shown  in  Fig.  2-5. 
For  example  an  ideal  relay  nonlinearity  can  be  seen  as  a  saturation  nonline¬ 
arity  which  linear  part  with  a  value  of  infinite  slope. 

The  operational  characteristic  of  a  nonlinear  element  which  operates  in 
a  control  system  depends  on  the  operating  condition,  just  as  a  vacuum  tube 
for  which  the  operating  characteristic  depends  on  the  operating  point  of  its 
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(a)  Saturation 


(b)  Variable  Gain 


(c)  bead.  £cne 


(d)  Hysteresis 


Fig.  2  5  3- sic  Concepts  of  Ilonlinoar 
Element 
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characteristic  curves.  As  a  saturation  nonlinearity  element  with  a  small 


forcing  function  input,  it  operates  linearly.  On  the  other  hand,  if  the 
forcing  function  is  very  large  and  with  a  high  frequency,  the  actual  opera¬ 
ting  characteristic  is  like  an  ideal  relay.  Hence,  the  actual  operating 
characteristic  of  a  nonlinear  element  may  be  said  to  be  a  function  of  the 
amplitude  and  frequency  of  the  forcing  function. 

2-6  Restoring  Force  Function  Investigation  of  a  Nonlinear  Element: 

Consider  the  block  diagram  of  Fig.  2-6,  in  which  the  output  is  a  func¬ 
tion  of  input,  usually  this  function  in  a  control  system  is  called  the 


Nonlinear 


e 


i 


e 


o 


Element 


Fig.  2-6  Basic  Relation  Between  Input 
and  Output  of  a  Nonlinearty 


m 


Restoring  Force  Function". 


Assume  the  input  is: 


eA  »  E(t) 


(2-23) 


The  output  will  be  a  function  of  input,  as: 


eQ  *  f(E)  (2- 

9 

From  the  numerical  analysis  and  curve  fitting  process,  the  function 
of  f(E)  can  be  expressed: 


(2-24) 


+  a  E 
n 


+  b  E 


n 


1/n 


(2-25) 
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Equation  (2-25)  can  be  in  the  general  form: 


f  (E)  «  a  +  v  •  a  Ei  +  b  El/i  (2-25a) 

°  t"«'  1  1  i  «^1 

for  this  paper,  the  nonlinear  characteristic  is  considered  symmetrical  to 
the  original  point,  that  means: 

f(-E)  «  -f (E)  (2-26) 

It  is  an  odd  function,  the  coefficients  of  even  order  of  equation  (2-25)  are 
zero,  it  becomes: 

3  5 

f  (E)  «  a^  +  a^E  +  a^E  +  a^E  + . * . . 

+  b^E  +  +  bjE1^5  + .  (2-27) 

For  the  saturation  and  variable  gain  nonlinearity,  J9  C.  West^  uses 
the  restoring  force  function  as  a  form: 

f (E)  *  axE  +  a3E3  (2-28) 

in  which  a^  is  greater  or  equal  to  zero,  and  a^  is  either  greater,  or  less, 
or  equal  to  zero  corresponding  to  the  case  of  a  "Hard  Spring",  "Soft  Spring", 
and  "Linear  Spring"  saturation  nonlinearity  characteristic  respectively. 

A  set  curves  for  all  values  of  a„  is  shown  in  Fig.  2-7. 


Fig.  2-7  Characteristic  curves  of  Eq.  (2-28) 
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In  order  to  generalize  for  the  condition  of  saturation  nonlinearity 
with  dead  zone,  equation  (2-28)  can  be  written  with  a  constant  as: 

f (E)  *  -aQSign  E  +  a^E  +  a^3  (2-29) 

if  the  dead  zone  is  less  than  one,  the  characteristic  curve  can  be  shown 
in  Fig.  2-8. 


r  u~> 

&3  rO 


Fig.  2-8  Characteristic  Curves  of  Eq.  (2-29). 

As  discussed  previously,  when  a  saturation  nonlinearity  is  operating 
with  a  larger  amplitude  than  the  saturation  voltage  and  a  high  frequency, 
the  operating  characteristic  is  like  an  ideal  relay.  On  the  other  hand, 
under  some  suitable  conditions  a  relay  nonlinearity  control  system  may  ex¬ 
hibit  some  sub-harmonic  forced  oscillations,  as  proved  by  the  paper  of  A.  M. 

12 

Hopkin  and  K.  Ogata  .  This  result  is  the  same  as  a  saturation  nonlinearity. 

Actually,  the  operation  of  a  relay  is  not  a  perfect  discontinuous 
characteristic  element,  when  the  relay  is  operating  from  static  to  closed, 
there  will  be  a  little  amount  of  time  delay  from  starting  to  close  con¬ 
tacts.  Therefore,  the  operating  characteristic  of  an  ideal  relay  may  be 
considered  a  saturation  nonlinearity  with  a  very  large  slope  of  the  linear 
part,  a  sketch  of  it  is  shown  in  Fig.  2-9. 
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Fig.  2-9  Operating  Characteristic  of 
an  Ideal  Relay 

From  the  above  discussion,  and  for  a  proper  value  of  constants,  the 
restoring  forcing  function  can  be  expressed  as: 

f  (E)  =  a^  +  a3E3  +  a^5  (2-30) 

Similarly,  if  we  consider  the  dead  zone  the  equation  of  (2-30)  can 
be  expressed: 

f (E)  =  -aQ  SignE  +  axE  +  a3E3  +  a^5  (2-31) 

13 

M.  J.  Abzug  investigates  the  restoring  forcing  function  of  an  ideal 
relay  nonlinearity  by  a  fifth  root  or  cube  root  of  the  input  as: 

f (E)  »  b5E1/5  (2-32) 

or; 

f(E)  -  b3E1/3  (2-32) 

A  more  detailed  discussion  about  frequency  response  by  using  the  re¬ 
storing  forcing  function  will  be  covered  in  the  next  two  chapters. 
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CHAPTER  III 


FUNDAMENTAL  FREQUENCY  FORCED  OSCILLATION  IN  NONLINEAR  CONTROL  SYSTEMS 
3-1  General  Description: 

Consider  a  system  shown  in  Fig.  3-1,  in  which  G^(s)  is  the  transfer 
function  of  a  controller  element.  *s  t*ie  transfer  function  of  a 

motor  and  gear  element,  G^  is  a  nonlinear  element.  H(s)  denotes  a  transfer 
function  of  a  feedback  network. 

When  a  periodic  forcing  function  is  applied  to  the  input,  the  output 
of  the  system  may  be  or  may  not  be  a  periodic  function.  If  it  is  a  peri¬ 
odic  output  and  the  frequency  is  the  same  as  the  frequency  of  the  input 
forcing  function,  it  is  said  that  the  system  is  operating  in  forced  oscil¬ 
lation  with  the  fundamental  frequency. 


V.  ! 

y 


^2. 


Fig.  3-1  Block  Diagram  of  a  Feedback 
Control  System 

There  are  some  methods  to  represent  the  forced  oscillation  in  a  non¬ 
linear  system.  If  we  keep  an  input  forcing  function  with  constant  frequency 
and  vary  the  amplitude  of  input,  a  response  curve  can  be  shown  in  the  " In¬ 
put  vs  Error"  plane,  Fig.  3-2,  in  which  there  is  a  cut  off  amplitude;  that 

means  when  the  input  is  smaller  than  0  ,  there  is  no  forced  oscillation. 

r-m 

The  other  representation  is  with  the  input  amplitude  constant  vary 
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Fig.  3-2  Forced  Oscillation  Represented 

in  the  E  vs  0  Plane 
r 

the  frequency  of  input.  It  is  shown  in  "Error  vs  Frequency11  plane.  Fig. 
3-2a  is  an  example,  it  is  actually  a  closed  loop  frequency  response  of  the 
system. 

tr 

\ 

/ 


\ 

Fig.  3-2a  Forced  Oscillation  Represented 
in  the  E  vs  ~0  Plane 

Either  represented  in  E  vs  0^  or  E  vs  co  plane,  both  of  them  can  be 
represented  in  a  phase  plane  (E  vs  E) .  A  typical  phase  plane  of  forced 
oscillation  from  the  computer  is  shown  in  Fig.  3-3,  in  which  it  is  represent¬ 
ed  by  an  ellipse.  The  size  of  it  is  changed  with  the  frequency. 

3-2  Basic  Equation  for  Forced  Oscillation  of  2nd  Order  Nonlinear  Feedback 
Control  System: 

Consider  a  2nd  order  feedback  control  system  shown  in  Fig.  3-4: 
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Fig*  3-4  Block  Diagram  of  a  2nd  Order 
Feedback  Control  System 

The  output  equation  in  Laplace  form: 


Vs> 


MID _ 

s(S+^  ) 


(3-1) 


in  the  differential  form: 

0  +  «.  9  =  Kf  (E) 

c  c 


(3-2) 


Where : 

ec(s)  *  er(s)  -  E(s)  (3-3) 

Equation  (3-2)  becomes: 

E  +  d  E  +  Kf (E)  =  0r  +  ^  0t  (3-4) 

Assume  the  input  forcing  function: 

0r(t)  »  F  Cos  (tt>  t  +  9)  (3-5) 


in  which  0 

is  the  phase  difference  between  the 

error  signal  and  forcing 

function  of 

input.  Inserting  equation  (3-5)  in 

to  equation  (3-4): 

E  +  E  +  Kf  (E)  *  -A  Cos  (it  t 

+  0)  — B  Sin(tO  t 

+  0)  (3-6) 

or: 

E  +  U-  E  +  Kf  (E)  -  C  Cos 

(i^  t+Jt  +  0  +  0) 

(3-6a) 

Where : 

> 

H 

6 

(3-7) 

B  »  oc  Os  F 

(3-8) 

„  / .  2  2 

C  *  A  +  B 

®  FtiJvC'  -r  4)-'- 

(3-9) 

0  *  Tan  *B/A  * 

Tan  ^  u  j 

(3-10) 
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Equation  (3-6)  or  (3-6a)  is  a  basic  differential  equation  of  a  2nd  order 
nonlinear  feedback  control  system,  in  which  f(E)  is  the  nonlinear  restor¬ 
ing  force  function,  and  K  is  a  constant  of  system. 

3-3  Forced  Oscillation;  Saturation  Nonlinearity  Without  Damping: 

Consider  a  system  without  damping,  i.e.^c  *  0,  equation  (3-6)  becomes: 

E  +  Kf (E)  «  -A  Cos(o>t  +  0)  (3-11) 

Let  the  restoring  force  function: 

f (E)  =  a^E  +  a3E3  (2-28) 


Equation  (3-11): 

E  +  K^E  +  a3E3)  -  -A  Cos(cu  t  +  0)  (3-12) 

7  8 

By  the  harmonic  iteration  method  *  ,  the  solution  of  equation  (3-12)  can 
be  expressed  as  a  series  of  odd  order  periodic  harmonics: 

e  (t)  *  Ex  Cos  ait  +  E3  Cos3^t  +  E$  Costa;  t  .....  (3-13) 

in  which  E^,  E^,  E^  ........are  the  amplitude  of  harmonics  of  e(t). 

If  we  consider  only  the  fundamental  frequency  forced  oscillation,  we 
can  assume  one  solution  is: 


e(t)  *  E^  CosuDt  (3-14) 

Where  E^  is  the  fundamental  frequency  amplitude  of  error  signal,  u)  is  the 
frequency  of  forcing  function  and  forced  oscillation;  is  to  be  determined. 
Since : 

E  *  -  oT  E  Cos  A)  t  (3-15) 

E  m  -  E  Sinu)t  (3-16) 

E3  *  ^  E^3(3  CosiOt  +  Cos3t'0t)  (3-17) 

Inserting  equation  (3-14)  into  (3-17)  in  equation  (3-12): 

2  3  3  1  3 

(-  00  Ex  +  Ka1E1  +  -  Ka3Ei  )  Cos  a)t  +  -  a^  Cos3^t  (3-18) 

*  -A  (Cos  aJ  t  Cos©  +  SinCjt  Sin0) 

neglecting  the  higher  order  harmonic  terms  and  equating  the  coefficients 
of  CosaJt  and  SinWt,  thus: 
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(3-19) 


(-»C2E,  +  Ka  ,E,  +  7  Ka  El'')  =  -  ACos© 

1  1  1  4  3 

ASlnut  Sln6  =0  (3-20) 

Equation  (3-20),  the  term  of  ASin  Ut  is  not  always  equal  to  zero,  the  term 

of  Sin0  should  be  zero.  Hence  the  value  of  6  will  be  either  zero  or  Jt. 

Inserting  equation  (3-7)  into  equation  (3-19): 

Ka  E  3  Ka  E  3 

- — - )  1/2  <3'21> 


For  0  =  0.  and: 


a)  = 


.Ka  El r 7 
_  (  1  4 


Ka3El‘ 


(3-22) 


for  6  *  *. 

As  a  check,  in  the  case  of  F  =  0,  i.e.  no  forcing  function  input,  equa¬ 
tions  (3-21)  and  (3-22)  give  the  free  oscillation  frequency  of  system: 

—  3  2  (6> 

l O  -  t'K  (al  +  jy  a3E  (3-23) 

For  the  linear  case,  i.e.,  a^  *  0,  the  natural  frequency  : 

u)ti  «  (ajK)'^  (3-24) 

It  is  exactly  the  same  as  we  discussed  in  the  linear  system. 

Recall  equations  (3-21)  and  (3-22);  the  value  of  a )  should  be  real, 
that  means: 

*1  +|-  a3E!2>  0  (3-25) 

for  both  the  phase  angle  is  either  zero  or  rt.  And  : 

Ej  >  F  (3-26) 

for  the  phase  angle  is  zero  only. 

Equations  (3-25)  and  (3-26)  are  the  conditions  for  existence  of  the 
fundamental  frequency  forced  oscillation. 

A  free  oscillation  frequency  response,  i.e.  F  *  0  can  be  sketched 
from  equation  (3-23),  any  value  of  force  function  other  than  zero  also  can 
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be  sketched  on  both  sides  of  the  response  curves  of  F  *  0.  It  is  to  be 
noted,  that  the  phase  angle  between  the  force  function  and  the  error  signal 
is  opposite,  i.e.,  0  *  n,  when  the  response  curves  are  to  the  left  of  curve 
for  F  =  0,  and  the  phase  angle  between  them  is  in  phase,  when  the  response 
curves  are  to  the  right  of  it.  In  other  words,  it  is  according  to  whether 
the  frequency  is  less  or  greater  than  the  frequency  of  free  oscillation  for 
that  particular  amplitude  of  error  signal  with  a  constant  of  force  function 
input . 

A  typical  error  vs  frequency  (E^ ^  cu  )  response  curves  and  phase  re¬ 
lation  sketch  for  a  different  value  of  forcing  function  and  different  char¬ 
acteristic  of  nonlinearity  are  indicated  schematically  in  Fig.  3-5.  The 
response  curve  for  free  oscillation,  corresponding  F^O,  is  drawn  as  a 
dash  line. 

In  this  respect  the  behavior  of  the  nonlinear  oscillation  is  the  same  as 
that  of  the  linear  oscillation.  One  sees  that  the  response  curve  in  the 
nonlinear  cases  could  be  thought  of  as  arising  from  those  for  the  linear  case 
by  bending  the  latter  to  the  right  for  a  hard  spring  saturation  nonlinear¬ 
ity,  and  to  the  left  for  a  soft  spring  saturation  nonlinearity. 

A  response  curve  in  " Input  vs  Frequency"  plane  (  ©c  ^  Cl)  )  can 
also  obtained  from  relation  of: 

ec(t)  =  er(t)  -  e(t)  (3-27) 

or: 

0  (t)  «  F  Cos(a)  t  +  0)  -  E,  Cos^t  (3-28) 

c  1 

In  the  case  of  0  =  0,  equation  (3-28)  becomes: 

0c(t)  =  (F  -  Ep  Cos  ujt  (3-29) 

Hence : 

0c  =  F  -  Ex  (3-30) 
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(o)  Hard  coring  S-iur  ti  n,  a-?  >  0 


© 

*  <jO 

0 


— i — 

<V /) 


(c)  Soft  spring  Sr.turatio  ,  <  0 

?  £" 

Fig.  Fespcnsc  Curves  for  the  S-.iure.tion 
x  I ‘online- rity  iithout  D..nping 
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or: 


E  *  F  -  9 
1  c 

Inserting  equation  (3-31)  into  equation  (3-21): 

K.  <6C  -F)  ♦  |  Ka3  (8C  -  F)3 

a }  *  ( - a - 

C 

For  the  case  of  ©  =  it,  equation  (3-28)  becomes: 

0  (t)  «  (F  +  E-)  Cos  <*)  t 
c  1 


Hence : 


1 


-  (f  +  ec) 


Inserting  equation  (3-34)  into  equation  (3-23): 


K  ,^8.  *  F)  -t  |  a3  (8C  f  F)3  m 
0 


(3-31) 

(3-32) 

(3-33) 

(3-34) 

(3-35) 


From  equations  (3-32  and  (3-35),  if  F  s  0,  the  result  is  the  same  form 
as  equation  (3-23),  that  means  the  response  curve  for  F  *  0  is  the  same  as 
plotted  in  the  CO  plane,  only  the  difference  is  that  the  value  of 

changes  to  0^.  The  response  curves  for  the  values  other  than  F  s  0  are 
changed,  and  the  phase  relationship  between  the  input  force  function  and  the 
output  is  reversed. 

There  is  a  numerical  example  for  plotting  the  response  curve  in  E  u) 
and  ©c  v  ^  plane  respectively  in  section  3-7  of  this  chapter,  and  also 
experimental  response  curves  are  plotted. 

From  the  analysis  of  the  above  and  the  numerical  example  in  section  3-7, 
it  is  shown  that  a  saturation  nonlinearity  feedback  control  system  can  have 
fundamental  frequency  oscillations.  This  not  only  depends  on  the  character¬ 
istic  of  the  nonlinearity,  but  also  depends  on  t he  amplitude  of  input  force 
function. 

3-4  Forced  Oscillations:  Saturation  Nonlinearity  with  Damping. 

As  explained  previously,  if  there  is  no  damping  present  in  the  system, 
there  is  an  oscillation  error  signal  e(t)  *  E^Cos(a)t)  either  in  phase  or 
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out  of  phase  of  180  degrees  with  the  forcing  function  of  inputs  In  the 
case  where  damping  is  present  however,  the  forced  oscillation  displacement 
and  the  forcing  input  can  be  expected  to  be  out  of  phase,  just  as  in  the 
corresponding  of  linear  system. 

Recall  the  basic  equation  (3-6)  and  the  restoring  force  function  equa 


tion  (2-28): 

E  +  'E  +  Kf (E)  *  -A  Cos(«'t  +  ©)  -  BS in (  ■  t  +■  ©) 
f (E)  *  a  E  +  a^E3 

(3-6) 

(2-28) 

Then  equation  (3-6)  becomes: 

E  +  v4E  +  Ka^E  +  Ka^E^  m  -ACos(^t  +  0^  -  BSin(  t  -  0) 

(3-36) 

Recall  one  solution  of  equation  (3-36)  is: 

e(t)  =  E  ^  Cos  u)  t 

(3-15) 

Inserting  equation  (3-15)  into  equation  (3-36): 

(-to2 E  +  Ka1E1  +  f  Ka  Ej3)Cos«,'t  -c^E^irWt  +  ~  Kb  E13Cos3  t 

»  (-ACos©  -  BSin©)  Cos^  t  +  (ASin©  -B  Cos©)  .iin  ■  t  (3-37) 
Neglect  the  higher  order  terms  of  harmonic  and  equates  the  coefficient 
of  Sin  u?  t  and  Cos^yt  respectively,  hence: 


(-CD2 E.  +  Ka  E  +  7  Ka  E, 3)  «  ACos©  -  BSin© 

1  11  4  3  1 

(3-38) 

E ^  *  A  SinG  -  B  Cos© 

(3-39) 

Squaring  equations  (3-38)  and  (3-39),  and  adding: 

(-  dtTE.  +  Ka,E  +  7  Ka,E  3)2  +  E  2  ■  A2  +  B2 

1  11431  1 

(3-40) 

or: 

3  32  2  2  2  2  f  2 

(-u,'El  +  Ka1E1  Ka^p  +  ^  cd  Ej  -  F  t-'  "  F 

(3-41) 

also: 

()f(E2  -  F2)  -  ^"[2E12(Ka1  +  ~  Ito  Ej2)  -U'  (E^  F2) 

3 

+  Ej,2  (Kax  +  4  ft  Ej2)2  =  0 

(3-42) 
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As  a  check  of  equations  (3-38),  (3-39)  and  (3-41);  first  consider 
the  response  curve  in  the  high  frequency  and  the  low  frequency  ranges, 
the  phase  relation  between  the  error  signal  or  control  signal  and  the  in¬ 
put  forcing  function  is  either  zero  or  180  degrees.  In  this  case,  the 
sine  terms  of  equation  (3-38)  and  (3-39)  are  zero,  the  response  curve  is 
the  same  as  equation  (3-20).  That  means  for  a  damped  system  if  the  phase 
angle  between  the  input  and  the  output  is  either  zero  or  180  degrees,  the 
frequency  response  is  exactly  the  same  as  the  system  with  no  damping. 

On  the  other  hand,  if  the  value  of  damping  is  very  small,  the  re¬ 
sponse  curve  will  be  very  close  to  the  sketch  of  those  of  Fig.  3-5.  The 
only  difference  when  very  small  damping  is  present  is  that  the  response 
curves  are  rounded  off  in  the  vicinity  of  curve  for  F  *  0.  A  sketch  of 
them  is  shown  in  Fig.  3-6. 

For  the  phase  investigation;  recall  equation  (3-39)  and  rewrite  in 
this  form: 


Sin0  -  JL  Cos©  *f>iE1/F 


(3-43) 


or : 


Sin(G  -  Tan"1--  )  = 


J-r  E  , 


F(  (3-44) 

Hence,  the  phase  angle  between  the  error  signal  and  the  forcing  function 
of  input:  ,  „ 

6  -  Sln‘l  !>v  i/2  +  Ta"'1  ir  <3-45> 

Examine  equation  (3-45)  in  the  case  of  *c  »  0,  ©  is  either  zero  or  jt, 
when  damping  is  present  the  value  of  ©  is  a  function  of  cO  and  if  the 
values  of  ^  and  F  are  fixed.  It  should  be  noted,  there  are  two  values  of 
phase  shift  for  one  value  of  error  signal;  one  is  for  the  low  frequency 
i)  and  the  other  is  for  the  high  frequency  one  is  from 

equation  (3-45)  and  the  other  is  from  the  180  degrees  minus  the  value  of 
equation  (3-45). 
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(a)  Linear  System,  a  =  0. 


(c)  Soft  Soring  Saturation,  a^<0. 

Fig.  3-6  Response  curves  for  The  S-tur-tir n 
Nonlinearity  ,-ith  .  a..ping 
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For  the  existence  conditions  for  forced  oscillation  when  damping  is 


present,  recall  the  equation  (3-42),  and  put: 


3  2 

Ka  4*  r  Ka0E-  »  X 
1  4  3  1 


and: 


(El2  -  F2)«j+  -[2E1ZX  *«;  (E^  -  F2)]'., 

Solve  equation  (3-47): 


2  2 

+  Ej  X‘  -  0 


(3-46) 


(3-47) 


U>  —  STtf.pjt  lZc‘A  (3-48) 

It  is  to  be  noted  that  the  conditions  for  forced  oscillations  in  an 
undamped  system  are  still  valid  for  a  damped  system  within  the  low  and 
high  frequencies,  in  other  words,  it  is  valid  for  the  condition  of  phase 
angle  either  zero  or  n. 

From  equation  (3-48)  where  forced  oscillation  exists*  the  value  of 
should  be  real,  hence;  the  condition  for  equation  (3-48)  to  be  real: 


and; 


or: 


2  2 
-  F  >  0 


2  2  2 
2E^X  -  ^(Ej  -  F  )  >  0 


X  ^ 


—Z—U  (E  2  -  F2)  >  0 
2E. 


(3-49) 


(3-50 


(3-51) 


and; 


or: 


[2E^X  -  ^  (E,2  -  r2)  J  2  -  4E1V(E1Z  -  F2)  >  0 


2.2,_  2  _2. 


4E12X2F2  +  (Ej2  -  F2)  >  ^Ej2X2  e*.  (Ej2  -  F2) 


(3-52) 


(3-53) 


Equation  of  (3-49);  (3-51)  and  (3-53)  are  the  conditions  for  the 
forced  oscillation  of  a  damped  system. 
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3-5  Jump  Resonance: 

Recall  a  sketch  of  frequency  response  curve  of  a  soft  saturation  non¬ 
linearity  shown  in  Fig.  3-7.  T7hen  a  experimental  frequency  response  is 
being  measured  for  a  system,  this  measurement  can  be  done  by  an  analog 
computer  to  measure  the  amplitude  of  the  input  and  the  output,  at  the  same 
time  measure  the  phase  relation  between  them. 


Fig.  3-7  The  Region  of  Jump  Resonance 
It  is  found  that  under  certain  conditions  of  amplitude  and  frequency 
an  infinitesmal  change  of  either  frequency  or  amplitude  of  the  input  signal 
causes  a  large  and  discontinuous  jump  in  amplitude  of  the  output,  at  the 
same  time,  a  discontinuous  phase  jump  occurs.  Fig.  3-8  shows  a  set  of 
amplitude  jump  and  phase  jump  of  a  typical  system  from  an  Electronic  Analog 
Computer  experimental  work. 

Fig.  3-7  is  shown  as  a  typical  system  with  a  soft  saturation  non¬ 
linearity  and  with  a  fixed  amplitude  of  the  input,  which  is  larger  than  the 
value  of  saturation  of  nonlinearity.  If  we  are  going  to  measure  the  frequency 
response,  being  with  a  low  frequency,  then  increase  the  frequency,  given  a 
response  from  A  to  B.  At  point  B  with  a  frequency  u)f  ,  an  infinitesimal 
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(b)  Lo'.rer  Jimp 


Fi;j  3  Q  A.  plitudc  and  Phase  Jimp 
of  a  Typical  Cjstera 
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increase  in  frequency  causes  the  amplitude  to  jump  upwards  from  point  B  to  C, 
and  with  a  phase  shift  from  below  90  degrees  to  a  value  above  90  degrees* 
Further  increases  in  the  frequency  cause  the  output  to  follow  the  gradual 
smooth  curve  to  D.  Now,  if  the  input  frequency  is  reduced  from  higher 
values,  the  output  will  follow  the  smooth  curve  DC.  When  the  frequency  is 
reduced  to  the  point  E,  if  the  frequency  is  infinitesimal  reduced,  the  out¬ 
put  jumps  down  from  point  E  to  point  F. 

It  is  found  that  either  from  the  theory  analysis  or  from  the  experi¬ 
mental  work,  a  hard  spring  saturation  nonlinearity  can  also  produce  the 
jump  resonance,  only  the  difference  between  them  is  the  jump  in  the  oppo¬ 
site  direction  from  each  other.  A  typical  system  for  a  hard  spring  satura¬ 
tion  nonlinearity  will  be  shown  in  the  experimental  section* 

From  the  above  discussion,  the  portion  of  the  curve  from  point  B  to 
point  E  cannot  be  obtained  experimentally,  usually  this  portion  of  the 
curve  corresponds  to  an  unstable  operating  condition* 

3-6  Lower  and  Upper  Jump  Frequency  of  the  Jump  Resonance: 

As  previously  discussed,  the  jump  frequency  of  is  called  "Upper 
Jump  Frequency'*  and  the  jump  frequency  of  6 OmZ.  is  called  **Lower  Jump  Fre« 
quency".  Recall  the  equation  (3-41)  and  the  sketch  of  the  response  curve 
shown  in  Fig.  3-7: 

(-^Ej  +  RBjEj  +  —•  Ka3E13)2  +  U*iV  E*  -  IT2  -#■ F2  (3-41) 

The  equation  of  the  locus  of  the  vertical  tangents  can  be  found  by 
differentiating  equation  (3-41)  implicitly  with  respect  to  and  setting 
d  /dE^  equal  to  zero,  the  result  is: 

(~^’“  +  Kfij  +  ~  Ka3E12)(-'o"+  Ka}  +  Ka  Ej2)  r^'-C  (3-54) 
Check  the  equation  of  (3-54),  if  the  damping  is  very  small,  the  term 
vC'uj  '  can  be  neglected,  that  leads  to  a  pair  of  equations: 

ur-  K(at  +  |  a^3)  -  0  (3-55) 
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a  K(aT  +  ~  a3E!2)  “  0  (3-56) 

The  curves  corresponding  to  these  equations  are  shown  in  Fig.  3-9. 

It  is  to  be  noted,  the  equation  of  (3-55)  is  the  same  as  the  response 
curve  equation  (3-23)  for  the  free  oscillation  of  an  undamped  system.  The 
equation  of  (3-56)  is  the  locus  of  the  vertical  tangents  of  curve  of  Fig. 
3-5. 

In  the  case  when  damping  is  present,  the  locus  of  vertical  tangents 
equation  (3-54)  will  be  directed  by  the  equations  of  (3-55)  and  (3-56), 
in  particular  there  should  be  one  branch  near  the  response  curve  for  free 
oscillation,  and  another  near  the  points  where  the  curves  for  the  undamped 
system.  In  other  words,  the  curve  must  appear  as  they  are  shown  in  Fig. 
3-10. 

Rearrange  the  equation  (3-54),  and  solve  for  id  ,  which  leads; 


zm,  J 'c*Wi  -T  3/ 


it  J7 


(3-57) 


or: 


0)  -  1  l  t  J  (-Lids  i 


3- 


58) 


In  which  .<,*  =  j,  Ka^,  it  is  the  natural  frequency  of  the  linear  system.  If 
we  know  the  constants  of  system,  and  the  amplitude  of  fundamental  fre¬ 
quency  forced  oscillation,  the  upper  and  the  lower  jump  frequency  can  be 
calculated  from  either  one  of  the  equations  (3-54),  (3-57),  and  (3-58). 


3-7  Analog  Computer  Analysis  for  a  Second  Order  System  with  Soft  Saturation 
Nonlinearity : 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  3-11. 
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FiC.3  9  Loci  oi  The  Vertical  Tangents 
ci  Response  Curve'  oC  A'-'.  Un.'.-^ipinr  h'onlinoar  Syctcn 


36 


t  J  - 


l 


(b)  ■  0 


,  .on: 


10  Looi  Oi  Vertical  lanocnv.s  of 
.ryes  of  a  I:\ir:ed  I'cniincar  Systen 
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i 

Fig.  3-11  Block  Diagram  of  a  2nd  Order  System 
with  a  Soft  Saturation  Nonlinearity 

In  which  the  restoring  force  function  of  nonlinearity  is: 

f (E)  -  2E  -  0.032s3  (3-59) 

The  characteristic  curve  of  restoring  force  function  is  drawn  in  Fig.  3-12 

where : 

E  *  1.00  1.50  2.00  2.50  3.00  3.50  4.00  5.00  6.00 

f(E)«  1.97  2.89  3.75  4.50  5.15  5.65  5.95  6.00  5.10 

Before  the  experimental  analysis  consider  the  theoretical  calculations 
and  curves*  First  consider  the  response  curve  on  the  B:  plane  for  the 

system  without  damping,  and  with  an  input  function: 

0^(t)  *  3  Cos(o;t  +  8)  (3-60) 

Recall  the  equations  of  (3-21);  (3-22)  and  (3-23)  and  leads: 

0)  -  (4  -  0.048Ej2)1/2  (3-61) 

for  free  oscillation  and: 


4E  -  0.048E  3  1/2 

"  -  <— — r  73 - >  <3-62> 

for  the  phase  angle  0  ~  :rt,  and: 


4E  -  0. 048E  2  1/2 

<-'t  -  (  - k - - -}  (3-63) 

E  -  3 

for  the  phase  angle  9=0. 
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The  calculations  and  curves  of  equation  (3-61)  to  (3-63)  are  shown 


in  Table  3-1  and  Fig.  3-13  respectively. 


Table  3-1: 

Calculation  for 

response  curves  on  the 

E ^  ''  plane  : 

Amplitude  of 

Frequency  for 

Frequency  for 

Frequency  for 

Error  Signal 

Free  Oscill. 

0  =  IT 

0^0 

(  £/) 

(  4  ) 

(  ) 

(  ) 

0.00 

2.00 

/ 

/ 

0.50 

1.996 

0.75 

/ 

1.00 

1.99 

0.99 

/ 

2.00 

1.95 

1.23 

/ 

3.00 

1.89 

1.37 

/ 

3.50 

1.84 

1.35 

4.88 

4.00 

1.80 

1.36 

3.60 

5.00 

1.68 

1.32 

2.65 

6.00 

1.53 

1.23 

2.16 

7.00 

1.29 

1.08 

1.73 

8.00 

1.00 

0.80 

1.22 

9.00 

0.35 

0.28 

0.42 

For  the  calculation  of  response  curves  and  phase  shift ,  when  damping  is 


present,  recall  equation  (3-42),  and  leads: 

2ei 

tltli.0.251-  — 2 — 2~  <4  ■  )]  + 

Ei  -3 


C1  2  2 
— tt — 9  (4  -0.048E,)  =  0 

E  ,2-  3 2  1 

(3-64) 


for  the  phase  shift,  recall  equation  (3-45) 

-1 


0  «  S  in 


Ei^ 


1/2  +  Tan 


1 


(3-45) 


F  (  <0^+  -L7-  )  & 

It  is  to  be  noted,  from  equation  (3-45),  and  the  response  curve  for  an  un¬ 
damped  system  shown  on  the  e,  -viO  plane,  when  the  frequency  is  larger  than 
the  frequency  of  free  oscillation  with  a  particular  amplitude  of  input, 
the  phase  angle  8  is  zero  and  is  rt,  when  the  frequency  is  below  the  free 
oscillation  frequency.  Hence;  the  phase  calculation  from  equation  (3-45) 
should  be  modified  as: 


0.25E, 


8  *  180 


(Sin 


-1 


0.25 


3(  *'+  0.252)1/2 


+  Tan 


(3-65) 
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for  i 


,  and : 


6  =  Sin 


-1 


0.25E 


1 


+  Tan 


-1 


0.25 


2  /  y  "> 

3 (  «  +  0.25  yc 


(3-66) 


for  tt  . 

If  the  response  curve  is  presented  on  the  plane,  the  condition 

of  equation  (3-65)  and  (3-66)  should  be  reversed. 

The  calculation  for  the  system  with  damping  of  0.25  and  the  curves  are 
shown  in  Table  3-2  and  Fig.  3-14  respectively. 

The  jump  frequencies  of  system  can  be  calculated  from  equation  (3-57) 
and  also  can  be  seen  from  the  response  curve.  It  is  shown  schematically 
in  the  diagram  of  Fig.  3-14. 


Table  3-2:  Calculation  for  Response  Curves  (with  a  damping  of  0.25) 


Amplitude  of 

Frequency 

Phase  for 

Frequency 

Phase  for 

Error  Signal 

for  46  <Lv'f) 

<  lOyi 

for^V't 

y  tsK/  -vj 

(  c-)  ) 

c  ^  > 

( 0  ) 

( to  > 

( i) ) 

1.00 

1.00 

164.00 

/ 

/ 

2.00 

1.24 

161.00 

/ 

/ 

3.00 

1.35 

158.50 

/ 

/ 

3.50 

1.36 

156.50 

4.83 

3.00 

4.00 

1.37 

155.00 

3.58 

7.00 

4.50 

1.40 

153.00 

2.80 

12.00 

5.00 

1.36 

151.50 

2.50 

15.50 

6.00 

1.24 

144.50 

2.05 

21.00 

7.00 

1.09 

135.00 

1.64 

28.50 

8.00 

0.87 

121.50 

1.16 

45.50 

8.70 

0.65 

90.00 

0.65 

90.00 

For  the  experimental  analysis, 

the  first 

investigating 

is 

the  restor- 

ing  forcing  function 

.  It  is  to  be  approximated 

as  a  two  slope 

lines,  one 

slope  is  2,  and  the 

other  is  zero. 

The  circuit 

setting  in  the 

analog  com- 

puter  is  shown 

below 

,  and  also 

all 

components 

to 

be  used  are 

schematically 

in  the  circuit 

of  Fig.  3-15. 

The  characteristic  curve 

for  the  restoring 

forcing  function  from  the 

computer  is  shown  in 

Fig.  3-16 

;  00 

is  shown 

the 

waveform  of 

input  (e^) 

and  output  (e  ) ,  (b)  is  shown  the  characteristic  of  the  output  to  the  input. 
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Fig*  3-15  Operational  Diagram  of  a  Function  ~ 

Generator  for  f(E)  «  2E  -  0*032E 

The  operational  diagram  of  analog  computer  setting  up  for  system  of 
block  diagram  in  Fig.  3-11  is  shown  in  Fig.  3-17,  in  which: 

e(s)  =  -  £  <S>J  0-67) 

ms)  =  _ L—  £  (S)  (3-68) 

r  / 

iO)  =  7^5“  63-  (s)  (3-69) 

From  equations  (3-67)  to  (3-69)  system  scaling  is  shown  below: 


where : 

o^S)  -  1.00;  oC?**  1.00; 

1.00; 

oCtX,  **  2.00 

wi 

as 

s 1  Rfl/R1  =  1.00 

ai 

=  1.00; 

Rfl  -  1.00; 

Rx  =  1.00 

W2 

= 

a2  Rfl^R2  =  1,00 

a2 

»  1.00; 

r2  =  1.00 

W3 

= 

Rf2C£1/a4  .  4.00 

a4 

=  0.50; 

Rj2  *  1.00 

W4 

m 

a3Ef2M4R3  *  4'00 

a3 

=  1.00; 

R3  *  0.50; 

W5 

= 

V  E5C£2  .2.00 

a5 

=  1.00; 

R5  *  0.50 

cf2  =  1.00 

W6 

7S 

a£R  /R,.  =  1.00 
or  o 

a6 

=  1.00; 

R,  -  1.00 

6 

Rf 3  s  ll0° 

in  which  R  in  megohm;  C  in  U  f 
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rig, 3 -/7  Operational  Diagran  of  Analog 
Couputer  setting  up  for  fasten 


The  results  fiom  the  analog  computer  experimental  work  Is  condensed 
in  Table  3-3;  and  the  response  curve  with  phase  shift  represented  on 
plane  is  shown  in  Fig,  3-18  respectively.  Also  the  jump  frequencies  are 
indicated  schematically  in  these  diagrams, 

A  theoretically  calculated  response  curve  in  a,  X  plane  for  comparing 
is  shown  with  the  experimental  response  curve  of  Fig,  3-18.  Comparing  the 
results  from  the  computer  experimental  work  with  the  results  from  theoretical 
calculation,  it  is  shown  that  the  results  are  very  close  within  high  fre¬ 
quency  and  low  frequency  range.  There  is  a  resonance  phenomena  between  the 
jump  frequencies. 
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Table  3-3  Data  for  Response  Curves  from  Analog  Computer  for  System 

,/  «  0.25;  K  *  2;  f(E)  *  2E  -  0.032E3:  6\(t)  -  3.0Cos(„t  t  9) 


(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude 

of 

Phase  < 

Tnput 

Error  Signal 

Error 

Output 

Outpui 

(0)  ) 

(  tS,) 

(  *  ) 

(  <£>  ) 

( <9 ) 

0.500 

0.38 

180.00 

3.20 

0.00 

0.628 

0.50 

180.00 

3.40 

0.00 

0.755 

0.75 

175.00 

3.60 

5.00 

0.880 

1.25 

168.00 

4.00 

9.00 

1.000 

1.50 

163.00 

4.50 

12.00 

1.130 

2.15 

158.50 

5.00 

24.00 

1.260 

3.00 

150.00 

5.40 

35.00 

1.260<+) 

9.00 

/ 

6.50 

/ 

1.380 

8.25 

25.00 

5.50 

140.00 

1.510 

7.30 

14.00 

4.50 

155.00 

1.640 

7.00 

5.00 

3.80 

159.00 

1.760 

6.50 

0.00 

3.40 

165.00 

1.890 

6.20 

0.00 

3.00 

172.00 

2.010 

6.00 

0.00 

2.60 

180.00 

2.140 

5.50 

0.00 

2.20 

180.00 

2.260 

5.20 

0.00 

2.00 

180.00 

2.390 

5.00 

0.00 

1.60 

180.00 

2.510 

4.80 

0.00 

1.40 

180.00 

2.830 

4.40 

0.00 

1.20 

180.00 

3.  140 

4.00 

0.00 

1.00 

180.00 

3.460 

3.80 

0.00 

0.80 

180.00 

3.770 

3.70 

0.00 

0.60 

180.00 

4.090 

3.60 

0.00 

0.50 

180.00 

(b)  Frequency  Decreasing: 

Frequency  of 

Amplitude  of 

Fhase  of 

Amplitude 

of 

Phase  o 

Input 

Error  Signal 

Error 

Output 

Output 

(to  ) 

( *?,) 

(  S'  ) 

(  0K  ) 

(  ) 

2 . 51 0 

4.75 

0.00 

1.50 

180.00 

2.200 

5.20 

0.00 

2.00 

180.00 

1.890 

6.20 

0.00 

3.00 

180.00 

1.760 

6.80 

0.00 

3.20 

173.00 

1.640 

7.  10 

0.00 

3.80 

168.00 

1.510 

7.58 

o.co 

4.50 

162.00 

1.380 

8.30 

12,00 

5.50 

156.00 

1.260 

9.20 

23.00 

6.80 

140.00 

1.130 

10.00 

34.00 

8.00 

135.00 

1.000 

11.50 

65.00 

9.50 

120.00 

0.880 

11.00 

/ 

8.70 

/ 

0. 880(~) 

0.75 

175.00 

4.00 

0.00 
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The  resonance  phenomena  in  Fig.  3-18  is  actually  similar  to  the  re¬ 
sonance  in  the  linear  system.  It  depends  on  the  value  of  damping  in 
the  system.  Another  experimental  result  for  the  same  system  and  same 
value  of  input,  only  changed  the  damping  from  0.25  to  0.5,  and  is  shown 
in  Table  3-4  and  Fig.  3-19  respectively.  In  this  case  the  resonance  pheno¬ 
mena  is  much  lower  than  the  case  of  a  damping  0.25. 

It  is  to  be  noted,  if  the  damping  of  a  system  is  increased  the  range 
of  jump  frequencies  is  decreased.  For  this  reason,  when  the  damping  in¬ 
creased  to  limit  value,  at  which  no  jump  phenomena  existed,  it  means  that 
the  system  becomes  a  linear  system.  On  the  other  hand,  if  the  system  is  an 
undamped  or  a  very  low  damping  system,  it  will  be  very  oscillatory  or  a 
limit  cycle  exists. 

The  response  curve  represented  on  "Output  vs  Frequency”  ) 

plane  can  be  calculated  from  equations  (3-32)  and  (3-35).  An  experimental 
response  curve  shown  on  the  £c »  &  plane  from  the  system  with  a  damping  of 
0.25  is  shown  in  Fig.  3-20. 

3-8  Analog  Computer  Analysis  for  a  Second  Order  System  with  Hard  Saturation 
Nonlinearity : 

Consider  a  system  with  a  block  diagram  shown  in  Fig,  3-21:  in  which 


A 


] 

.) 

A 

*  s 

\ 

Fig.  3-21  Block  diagram  of  a  Second  Order  System 
with  a  Hard  Saturation  Nonlinearity 

assume  the  restoring  forcing  function  of  the  nonlinearity  is: 

f (E)  =  0.425E  +  0.005E3  (3-70) 
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Table  3-4,  Data  for  Response  Curve  from  Analog  Computer  for  System 

^  =  0.50;  K  =  2;  f(E)  -  2E-0.032E3;  6,(t)  =  3. OCos (a) t+o) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

( *  > 

( > 

( *  > 

( *<.) 

(  C  ) 

0.628 

0.50 

173.00 

3.40 

7.30 

0.755 

0.75 

165.00 

3.50 

12.00 

0.880 

1.00 

155.00 

4.00 

16.50 

1.000 

1.50 

140.00 

4.75 

24.00 

1.130 

1.75 

125.00 

6.00 

30.00 

1.260 

8.30 

91.00 

8.00 

85.00 

1.380 

8.10 

47.00 

6.25 

120.00 

1.510 

7.50 

36.00 

5.25 

135.00 

1.640 

7.00 

25.00 

4.25 

148.00 

1.760 

6.50 

20.00 

3.50 

160.00 

1.890 

6.00 

12.00 

2.75 

172.00 

2.010 

5.50 

5.00 

/ 

/ 

2.140 

5.00 

0.00 

/ 

/ 

2.260 

4.50 

0.00  . 

1.75 

180.00 

2.390 

4.25 

0.00 

/ 

/ 

2.510 

4.00 

0.00 

1.25 

180.00 

2.810 

3.75 

0.00 

1.00 

180.00 

3.140 

3.50 

0.00 

0.75 

180.00 

3.460 

3.35 

0.00 

0.50 

180.00 

3.770 

3.25 

0.00 

0.30 

180.00 

4.090 

3.10 

0.00 

/ 

/ 

(b)  Frequency  Decreasing 

Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input: 

Error  Signal 

Error 

Output 

Output 

(  «>  ) 

(  e--\  ) 

(  0  ) 

(  O,  ) 

(  6  ) 

2.510 

4.30 

0.00 

1.25 

180.00 

2.200 

5.00 

0.00 

1.75 

180.00 

1.890 

5.75 

0.00 

2.75 

172.00 

1.760 

6.25 

7.50 

3.50 

160.00 

1.640 

7.00 

15.00 

4. 10 

148.00 

1.510 

7.75 

27.00 

5.00 

132.00 

1.380 

8.25 

32.00 

6.50 

120.00 

1.260 

8.50 

60.00 

7.50 

97.00 

1. 130 

7.00 

80.00 

8.20 

40.00 

1.000 

2.00 

130.00 

4.75 

19.50 

0.880 

1.25 

145.00 

4.00 

0.00 
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Fig.  3 -SI'1  Characteristic  Curve  of 
torin^  function :  f  )  —  0  J!-S5f'v0 . 003^3 


The  characteristic  curve  of  equation  (3-70)  Is  shown  in  Fig*  3-2la,  where: 

E  «  1.00;  2.00;  3.00;  4.00;  5.00;  6.00;  7.00;  8.00. 

f (E)  =  0.43;  0.89;  1.42;  2.02;  2.75;  3.63;  4.69;  5.96. 

A  theoretical  analysis  can  be  calculated  from  equations  (3-70)  and 
(3-21)  to  (3-23),  as  we  did  in  section  3-7.  Only  the  analog  computer  anal¬ 
ysis  is  considered  in  this  section. 

For  the  simulation  of  the  restoring  forcing  function  of  nonlinearity 

* 

from  the  analog  computer,  it  is  approximated  by  two  straight  lines,  one 
with  a  slope  of  0.425,  and  the  other  with  a  slope  of  2.00.  The  circuit  set¬ 
ting  in  the  computer  is  shown  in  Fig.  3-22,  and  all  components  to  be  used 
are  schematically  in  the  diagram. 


/  c 


r,  ^  ' 

-^1  .i  h  i  I- 


/Z4  ! 
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Fig.  3-22  Operational  diagram  of  the  Function  Generator 
for  f (E)  »  0.425E  +  0.005E3 

The  characteristic  curve  of  restoring  forcing  function  from  the  com¬ 
puter  is  shown  in  Fig.  3-23. 

The  operational  diagram  of  analog  computer  set  up  for  the  whole 
system  is  mostly  the  same  with  the  diagram  of  Fig.  3-17,  the  only  dif¬ 
ference  is  the  function  generator  circuit. 
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The  results  from  the  analog  computer  are  condensed  in  Table  3-5, 
and  the  response  curves  with  phase  shift  shown  on  the  rpt  and  S‘,  v  A 
plane  are  shown  in  Fig.  3-24  and  3-25  respectively.  In  which  the  input 
forcing  function  is  8  (t)  =3.50  Cos(^t  +  8). 
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Table  3-5.  Data  for  Respoi  se  Curves  from  Computer  for  a  System  with  a 
Hard  Saturation  Nonlinearity,  in  which:  =  0.25;  K  «  2; 

f(E)  =  0.425E  4  0.005E3;  9r(t)  =  3.5  Cos(^  t  4  0). 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

( to) 

(  £>  ) 

(  0  ) 

(  <S>c) 

(  <p  ) 

0.500 

1.75 

180.00 

5.50 

0.00 

0.628 

3.50 

180.00 

6.50 

4.00 

0.754 

6.50 

175.00 

9.00 

7.50 

0.880 

8.00 

168.00 

10.50 

10.00 

1.000 

9.00 

164.00 

12.00 

11.50 

1.130 

10.50 

158.50 

13.00 

13.00 

1.260 

12.25 

152.0Q 

14.50 

15.00 

1.380 

14.00 

145.50 

16.30 

17.00 

1.510 

16.25 

141.00 

18.75 

20.00 

1.640 

19.00 

134.00 

21. 30 

24.00 

1.760 

22.25 

124.00 

22.00 

37.00 

1.890 

26.00 

102.00 

25.00 

45.00 

1.950 

26.50 

/ 

25.50 

55.00 

2.010 

26.00 

24.00 

25.00 

90.00 

2.140 

4.50 

20.00 

0.86 

150.00 

2.260 

4.35 

15.00 

0.65 

165.00 

2.510 

4.25 

7.50 

0.50 

180.00 

2.830 

4.15 

0.00 

0.40 

180.00 

3. 140 

4.00 

0.00 

0. 10 

180  no 

^.770 

3.80 

0.00 

/ 

/ 

(b)  Frequency  Decreasing: 

Frequency  of 

Amplitude  of 

Phase  of 

Amplitude  of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(<D  ) 

(  £,) 

(  <P  ) 

(  £><  ) 

(  A  ) 

2.510 

4.25 

5.00 

0.50 

180.00 

2.400 

4.35 

10.5'> 

0.63 

180.00 

2.270 

4.50 

13.50 

0.68 

172.50 

2.140 

4.60 

15.00 

0.75 

170.00 

1.890 

4.85 

21.00 

1.13 

157.50 

1.760 

5.00 

23.50 

1.25 

155.00 

1.640 

5.35 

25.00 

1.60 

152.00 

1.510 

5.75 

26.00 

2.  13 

150.00 

1.380 

6.75 

27.50 

16.30 

18.50 

1.320 

12.25 

/ 

/ 

/ 

1.260 

11.70 

150.00 

14.50 

15.50 

1.130 

10.25 

155.00 

13.00 

13.50 

1.000 

9.00 

160.00 

11.70 

12.00 

0.880 

7.50 

172.00 

10.50 

10.00 
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CHA1TER  IV 


SUB-HARMONIC  OSCILLATION  IN  NONLINEAR  FEEDBACK  CONTROL  SYSTEM 
4-1  General  Description 

Up  to  the  preceding  chapter,  there  is  only  discussed  the  forced 
oscillations  for  which  the  frequency  is  the  same  as  that  of  the  external 
force,  A  nonlinear  control  system,  under  suitable  conditions  may  have  a 
steady  state  oscillation  in  which  the  main  component  has  a  frequency  which 
is  a  fraction  1/n  (n  is  a  real  and  positive  integer)  of  the  forcing  func¬ 
tion  frequency.  These  oscillations  are  called  MSub-Harmonic  Oscillation”. 

Recall  a  linear  system  (Section  2-3),  there  are  two  oscillation  terms 

in  the  transient  response,  if  the  frequency  of  the  free  oscillation  is 

C/)  /n  (n  is  an  integer),  then  a  force  function  of  frequency  <'/  can  excite 

the  free  oscillation  in  addition  to  the  forced  oscillation  of  frequency  c  . 

~/t/2 

But  as  t  increases,  e  '  damps  the  free  oscillation  until  only  the 
steady  state  oscillation  remains.  Hence  there  is  no  sub-harmonic  or  super¬ 
harmonic  oscillation  existing  in  a  linear  control  system  at  steady  state. 

In  the  case  of  nonlinear  control  system,  the  steady  state  response 
may  not  have  the  same  frequency  of  oscillation  as  the  frequency  of  the 
forcing  function.  In  fact,  the  frequency  of  the  system  response  depends 
not  only  upon  the  frequency  of  forcing  function  but  also  upon  the  amplitude 
of  the  input.  The  principle  of  superposition  no  longer  holds  and  the  linear 
circuit  characteristic  is  not  applicable.  In  general,  it  can  be  said  that 
the  nonlinear  system  can  have  wide  variety  of  almost  periodic  oscillations, 
the  frequencies  of  which  differ  from  the  frequency  of  the  applied  force  and 
vary  with  time  as  well  as  with  different  initial  conditions. 

Sub-harmonic  oscillation  always  existed  with  a  frequency  higher  than 
the  frequency  of  free  oscillation  with  a  particular  amplitude  of  error  signal. 
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The  response  curve  which  is  a  part  of  an  ellipse  or  a  hyperbola  in  the 

plane  is  always  nearly  parallel  to  the  right  side  of  the  free  oscil¬ 
lation  response  curve.  A  typical  sketch  of  a  soft  saturation  nonlinearity 
is  shown  in  Fig.  4-1. 


Fig.  4-1  Fundamental  and  Sub-harmonic  Response 

Curves  for  a  Soft  Saturation  Nonlinearity 

Under  some  suitable  conditions,  there  may  exist  any  order  of  sub¬ 
harmonic  oscillation,  or  only  one  type  of  order  can  exist  depending  on 
the  characteristic  of  the  nonlinearity  and  the  input  forcing  function. 

Sub  harmonic  oscillation  and  response  of  nonlinear  vibratory  systems 
are  often  discussed  in  the  field  of  mechanical  system  9  There  is  no 
body  of  knowledge  applied  to  feedback  control  servome chan ism.  For  the  feed¬ 
back  control  system,  there  is  a  standard  differential  equation  (3-6a),  it 
is  very  similar  to  the  equation  for  the  mechanical  system,  the  only  dif¬ 
ference  is  that  the  forcing  function  of  a  feedback  control  system  is  a  func¬ 
tion  of  input  frequency. 

Although  any  order  of  sub-harmonic  oscillation  can  exist  in  a  non¬ 
linear  feedback  control  system,  in  this  paper  only  the  sub-harmonic  of  order 
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2  and  3  are  discussed. 


4-2  Sub-harmonic  oscillation  of  Order  2  without  Damping: 

Sub-harmonic  oscillation  of  order  2  can  exist  in  most  nonlinear  feed¬ 
back  controls.  When  it  occurs,  the  frequency  of  output  (Co)  or  error 
(e)  is  only  1/2  the  frequency  of  input  (  6Jr  ) .  A  typical  diagram  from 
the  computer  experimental  work  shown  on  the  "Input  vs  Output"  plane  and  the 
waveforms  of  "Input  and  Output"  are  shown  in  Fig.  4-2  and  4-3,  respectively. 

Recall  the  basic  equation  for  forced  oscillation  in  nonlinear  system 
with  un-damping  Eq.  (3-11). 

E  +  Kf (E)  =  -A  Cos (  t  +  0)  (3-11) 

For  the  same  type  of  restoring  force  function  with  chapter  3,  leads: 

E  +  Ka1E  +  Ka^E3  =  -  A  Cos(<u  t  +  0)  (3-12) 

Assume  one  solution  of  equation  (3-12)  is: 

E  =  E  ^Cos  +  Ej  Cos  (*j  t  (4-1) 

Where  E^0  is  the  amplitude  of  sub-harmonic  of  order  2,  E1  is  the  amplitude 
of  fundamental  frequency,  since: 


E  =  4 


'1/2 


Cos  yiOt  - 


iV  Ej  Cos  u)  t 


(4-2) 


e4*  i-%4.  r  r 


,  .  }  *  1  3 

r  4  r  j  ly_  r  1  u, '  2  '~c  T  jf 


'■  K'<  r  i  -  '  r  -■  c<-  r  (4-3) 

We  are  only  interested  to  the  terms  of  second  sub-harmonic  and  fund¬ 
amental  frequency,  inserting  equations  (4-1)  to  (4-3)  into  equation  (3-12) 
the  result  is: 

(-  §‘:6yv  t  t'A.tfr,  r/F  1 1  fij,  G?  )  a.  J;+t  t  (~  u>"£,  i 


„  I 

r  1-,  r/  Y  L  J  '"o- 


-  ,wr 


A’  ' 


(4-4) 
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K  •"  4;  ^-0,2 5 

er(t)^'2.50Cos(  m.  t  -  o) 
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As  discussed  in  Chapter  3,  Section  3-3,  if  the  frequency  is  larger 
then  the  frequency  of  free  oscillation,  the  phase  angle  G  is  zero  on  the 
"Error  vs  Frequency"  plane,  hence  equation  (4-4)  leads: 

uor  —  4  -A-',  t  ,.fiy  *  L  (4*5) 


CO' 


r  A  l  Cz y  t  j  /-y 

E-  F 


(4-6) 


Equation  (4-5)  is  the  condition  for  existence  of  the  sub-harmonic  of  order 


2.  Solve  Ej/2: 


‘  L  ( 


-  ±.  Ik  T  6B^yyt 

=>  c/3  L  u>  s 

Since  E.^  stlou^  real;  which  leads: 


(4-7) 


00  <  (4-Kd,  *  ! &  Kas£fJ*'  a,  Z  v 

For  the  first  iteration  method,  consider  equal  to  zero: 

Equations  (4-5)  and  (4-6)  leads: 

O0~  -  Ah*-’-! 


(4-8) 


(4-9) 


E,  -  ft  (4-10) 

Inserting  equations  (4-9)  and  (4-10)  into  equation  (4-8): 

cjj  ^  t  &  if  a 3  <  o  (4-ii) 

From  Equation  ('-5),  it  represents  an  ellipse  or  a  phperbola  in  the  %y/2 
vs  o>  plane  depending  on  the  sign  of  a  also  to  has  a  maximum  when  a. 
is  less  than  zero,  and  a  minimum  when  a^  is  greater  than  zero.  Actually, 
the  sub-harmonic  oscillation  only  existed  within  an  interval  of  frequency. 
That  means  only  one  part  of  ellipse  or  hyperbola  can  exist,  as  shown  in  Fig, 


(4-4). 

4-3  Sub-harmonic  Oscillation  of  Order  2  with  Damping: 

Recall  the  basic  equation  (3-41)  for  a  damped  system: 

-  3 

E  -t  E  +  Ka^E  +  Ka^E  =  -  ACos(o;  t  4*  G)  -  BSin(at  +  0)  (3-41) 
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Assume  one  solution  of  equation  (3-41)  is: 


t-  —  Cr 


•  i  E,n  c-  <■  c,f.  !•/--?  a  t 


(4-12) 


Since : 


£  -•  -  J-AV.  1  /./y/  ‘i-v.x&t  r  zr,^ 

23  ’  jj:  o:  ' b, ^ /-tJ'V^/Afr  ■  £,fi S. >/ 2^6 

i ^  +■  i.  /-  //j  V-2.  *  i'-'/i'-fyi )  (U" 

t  ,  r  9rb,li  <~  .t-'/L  E,/<  )  (:<  5  ^ ' 

'  (^~'-iS  r  Tt7//>  r  ^ 


(4-13) 

(4-14) 


(4-15) 


Inserting  equations  (4-12)  to  (4-15)  into  equation  (2-41)  and  equal¬ 
ing  the  coefficients  of  Cos  *  Cos:C£,  and  ^ink^r,  give  a  result  of 

+J;fy7r.bW')  ““  0  (4-16) 

~£//ib~  LV  t-KH,  fzM  rt-Cl>(.fZ,A  r  r  ) 

^  Cot  o  [4  ■jsi')  $  (4-17) 

-k,\jhr  -  in,/,  *j  t  Ktit  9,i)  r  A'? 5  f-j,  r.y  rk£/J-ttb> 

-  Q.  ~  Q  C  (>  V  (P  (4- 16) 


From  equation  (4-18),  it  is  the  coefficient  of  the  term  of  Cos  ^  4,  t. 
where  Ei/2  ^^^ence' 
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(1)  Func.auentsl  Frequency. 


(a)  Soft  Saturation ;  a.j  -  0 


(B)  Ifcir'd  C^-BUi’-TiiC  ii  |  >  0 


Fig.  4-4  Response  Curves  for 
Sub  hnmcnic  Cscillatior.s 
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2 


(4-19) 


Solve  for 

*1/2  .  -  tv,  <  ■  <*-»> 

The  value  of  ^^{2  s^ou^  rea^  therefore,  the  existence  of  2nd  order 
sub-harmonic  with  damping: 

^  V  -Ml,  +  l  ^  ^  ^//i  f  U^i..  a3  0  (4*21) 

From  equations  (4-16)  to  (4-18)  for  the  first  iteration  method,  for 
a^  equal  to  zero:  hence: 

£u  '  *  4Ka ^  (4-22) 

(Ka  ^  ) E 1  +  a  •  Ej^  *  -  FCosG  -  FSinG  (4-23) 

(KaL  -  i)S  )ElB  -  AO;  Eu  -  to  FSinG  -  FCosG  (4-24) 

Solve  Ej^  and  from  equations  of  (4-22)  to  (4-24): 


E 


1A 


~M( 


. » i  j  i 

Z.  i  +-  * 


'IB 


A  /-7 


j,  ,  i. 

/ 


(4-25) 

(4-26) 


Inserting  equations  (4-25)  and  (4-26)  into  equation  (4-21).  A 
general  equation  for  existence  of  2nd  order  sub-harmonic  oscillation 
with  damping: 


03  =2.  <2  '!■  f-  Kab  h*L  O  -  pifr-t'*  C  i~  ^  ^ 

T  ~ 


/w,t  u ~  fn,  x  *  4 


V 
,  ; 


■-  f(^l*  .  ^  Cfr*  I  S1  A  j  f 


n* 


t  4  J 


(4-27) 


For  the  case  of  G  equal  to  0  or  «,  i.e.,  SinG  equal  to  zero  equation 
(4-27)  becomes: 


.V 


2- 


•"•yV  • 


r  A^;i 


(4-28) 


69 


From  equations  (4-11)  and  (4-28)  they  were  shown  that  either  the 
system  with  damping  or  without  it  the  second  order  sub-harmonic  oscilla¬ 
tion  only  existed  with  a  frequency  less  than  the  frequency  of  two  times 
the  natural  frequency  of  linear  system  for  the  case  of  less  than  zero, 
and  greater  than  this  frequency  for  the  case  of  a^  greater  than  zero. 

4-4  Sub-ha rmonic  Oscillations  of  Order  3  without  Damping: 

Typical  waveforms  with  jump  phenomena  of  3rd  order  sub-harmonic 
forced  oscillation  is  shown  in  Fig.  4-5;  (a)  shows  the  transient  from 
the  2nd  order  to  the  3rd  order  sub-harmonic;  (b)  is  shown  reversing. 

For  investigating  the  existence  conditions  of  3rd  sub-harmonic  forc¬ 
ed  oscillation  recall  the  basic  equation  of  (3-11);  and  leads: 


E  +  Ka jE  +  Ka3E  *  -  A  Cos(wt  +  G)  (3-12) 

Assume  one  solution  of  equation  (3-12)  is: 

E  ®  E1/3  Cos  ^  *,  t  +  Ej  Cost  t  (4-29) 

in  which  is  the  amplitude  of  3rd  order  sub-harmonic  oscillation^ 

hence: 

E  *  -  ^  El/3  ^os  *3  ^  ~  a  El  Cos  ^ t  (4-30) 

3  3 

E  =  *  EL/3<e21/3  +  E1/3E1  +  2EP  Cos3  J  C 


13  2 

"j  (E  ^  j  ^  ^  j  ^  E  ^  +  3E  ^  )  Cos  (v  t 

From  equation  (4-29)  to  (4-31)  and  equation  (3-12)  leads: 


(4-31) 


CK&-  j  '  f  'u' 


(  '7-' 


(4-32) 


-  *  •'  rjf.  *■'  •  •  r  ' 

Equation  (4-32)  E^  Is  not  equal  to  zero,  hence: 


>  i  a 


Z-V  t  tr  tt, 

/  7 


t 


(4-33) 


(4-34) 


Solve  from  equation  (3-34), 
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i*  is.  4-5 


‘..aveforns  of  Upper  and  Loner  Jimp  Transient 
for  A  Topical  3~d  Order  Sub-harnonic 
"  o  reed  0  sci  llation 
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A 


(4-35) 


The  value  of  must  be  real,  then: 

£  W.  r  «  7£,  '  -?  t  (4-36) 

Equation  (4-36)  is  the  condition  for  the  existenceof  sub-harmonic 
oscillation  of  order  1  of  a  without  damping  system: 

From  equations  (4-32)  and  (4-33),  by  the  first  iteration  method, 


for  «  0;  then: 
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9  Ka . 


(Ka^  -  6c  )Ej  *  F 

Solve  E^  from  equation  (4-37)  and  (4-38): 

E  s  -  F 

hl  8  b 

Inserting  equation  (4-39)  in  equation  (4-36)  the  condition  for  a  3rd 
sub-harmonic  oscillation  becomes: 


(4-37) 

(4-38) 

(4-39) 


U) 


?  •/’  -v  "  '  ~  x  it  K(iif ' 


Also  the  response  equation  of  equation  ^-3/  becomes: 


(4-40) 


(4-41) 


The  equation  of  (4-41)  represents  an  ellipse  or  a  hyperbola  in 
the  E^^  vs  ^  plane  depending  on  the  sign  of  also  has  a  maximum 

when  0,  and  a  minimum  when  >  0,  A  sketch  of  a  3rd  sub¬ 

harmonic  forced  oscillation  is  shown  in  Fig,  4-4, 

4-5  Sub-harmonic  Oscillation  of  Order  3  with  Damping: 

Recall  equation  (3-41): 

E  +  E  +  Ka^  +  Ka3E3  =  -  ACos(.4t  +  0)  -  3Sin(iv  t  f  0)  (3-41) 

Assume  one  solution  of  equation  (3-41)  is: 


E  *  E ,  1 ~Cos  *0  t  4*  E.  .Cos  ^  t  *4  E«_Sinw/  t 
1/3  3  1A  IB 


(4-42) 
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Since : 


«  1 

E  *  E.  ,  Q  Sin  —  *  t  -  E-  .  ^  Sin  t  +  E.  0  Cos  t 

1/5—3  1A  IB 


6  -  '  |e1/3  cos  |  t 


E.  .  .v  Cos  t  -  -  E. _  Sin  t 
1A  IB 


(4-4  .) 
(4-44) 


'3  -  I  El/3(El/32  +  Sl/3E1A  *  2EU2  *  2ElB2)COS  3  * 


4El/3  ElBSin3  U  t  +  1/,4^El/3  +  El/3  E1A  +  3E1A  + 

3eiaeib'>  CosWC  ♦  7  E1B  «|»2  *  2  El/32  *  E1A2>  **»“  *  ••<*-«> 

4 

From  equation  (4-42)  to  (4-44)  and  (3-41) ,  by  equaling  the  coefficients 
of  Cos  %  t,  Sin  t,  Cost  t,  leads: 


Mi,-  r; 


^  1 


2  c  t  r  > 


(4-46) 


/  t  r  . 


■  .  O',,  !  »  )  tc/^  ,  - 
—  w  -  4  -  f- 


^  A  ~  r  *  . 


'  "fA'  -  <5?  7 "■*>: 


h 


(4-47) 

(4-48) 


Rearrange  equation  (4-46)  and  solve  it  for  *-n  terms  °£ 


eia*  an,i  eib: 

% 


'  7  ,  ^ 


//'/  > 


'/r* 


* 


here  ^  should  be  real*  hence: 


7  '  ^ 

or 


/ 


r  ;> 


or: 


(4-49) 


(4-30) 


,  .  3  A‘<:  4  / 


r  <<  <  " 

!  *?/ 


(4-51) 


Equation  (4-51)  is  the  condition  of  3rd  sub-harmonic  oscillation 
existing  in  saturation  nonlinearity  with  damping  system. 

By  the  iteration  method,  equation  (4-46)  to  (4-48)  leads: 

J  «  9Kax 

(Ka^  -  cv  )E^  +  U.  X  »  -  a  FCos  8  -  ^  FSin© 


(Ka,  -  A  )E._  -  v  E. .  =  i  FSin  0  -  y,  FCos© 
N  l  IB  1A 


(4-53) 

(4-54) 

(4-55) 
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Solve  equation  (4-53)  to  (4-55), 


E  =  ~l 
1A  J 


(4-56) 


'IB 


(4-57) 

L  f  * 

Substituting  equations  (4-56)  (4-57)  into  equation  (4-52)  and  be¬ 
comes  in  a  general  form: 


i)  Z 


,  L  L  t  * 


L*  >*t'V  ^ 


-  f  t 

i  <  r 


e ■  !  i  L3s£  i 


-r  -v  r 


(4-58) 


In  the  case  of  Q  equal  to  0  or  n,  i.e,  Sin8  equal  to  zero  equation  (4-58) 
becomes : 


a  i  3 


lJ  ^f\ 


A 


t* 


/  V 


/  2- 


y\ 


(4-59) 


In  conclusion,  the  condition  for  existence  of  sub-harmonic  oscilla¬ 
tion  in  a  nonlinear  feedback  control  system  depends  on  the  characterics 
of  the  nonlinear  element  and  the  amplitude  of  input •  The  response  curve 
of  any  order  sub-harmonic  oscillation  with  a  saturation  nonlinearity  is 
either  an  ellipse  or  a  hyperbola  which  depends  on  the  characteristic  of 
nonlinearity*  Also  the  frequency  of  a  sub-harmonic  oscillation  existed 
may  be  either  a  maximum  when  a  soft  saturation  nonlinearity  is  presented 
or  a  minimum  when  a  hard  saturation  nonlinearity  is  presented*  The  limit 
value  of  minimum  or  maximum  is  almost  equal  to  the  number  of  sub-harmonic 
order  times  the  natural  frequency  of  linear  system. 

4-6  Analog  Computer  Analysis  for  a  2nd  Order  Sub-harmonic  Oscillation 
System: 

Consider  a  system  with  a  block  diagram  shown  in  Fig.  4-6 1 
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Fig.  4-6  Block  Diagram  of  a  2nd  Sub-harmonic 
Oscillation  System. 

Assume  the  restoring  force  function  of  nonlinearity  is: 

f (E)  =  2E  -  0.0048E3  (4-59) 

The  characteristic  curve  for  equation  (4-59)  is  shown  in  Fig.  4-7,  in 
which: 

E  *  1.00;  2.50;  5.00;  7.50;  10.00;  12.50;  15,00. 

f (E)  =  1.995;  4.295;  9.400;  12.97;  15.20;  15.40;  14.80. 

For  the  simulation  of  restoring  force  function  of  nonlinearity  from 
computer,  it  is  approximated  from  two  straight  lines,  one  of  which  with 
a  slope  of  2;  and  the  other  with  a  slope  of  zero.  The  circuit  setup  in  the 
computer  is  identical  with  the  circuit  shown  in  Fig.  3-15.  The  only  dif¬ 
ference  is  the  value  of  bias  voltage  to  be  used.  The  characteristic  curve 
of  restoring  force  function  from  the  computer  is  shown  in  Fig.  4-8. 

An  operational  diagram  setup  in  the  analog  computer  is  almost  the 
same  with  the  diagram  of  Fig.  3-20,  only  one  difference  is  that  the  value 
of  constant  to  be  used. 

The  results  from  the  computer  are  condensed  in  Table  4-1  and  Fig.  4-9. 
Only  the  response  curve  is  of  interest,  hence,  only  one  response  curve  in 
the  error  vs  frequency  plane  is  presented.  It  is  shown  that  there  are  two 
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rig. 4- 7  Characteristic  Curve  Plotted 
of:  f(E)— 2E  -  0.00'!-oE3 
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jumps,  one  is  the  fundaments  1  frequency  jump,  and  the  other  is  the  second 
order  sub-harmonic  frequency  jump.  These  jump  phenomenas  are  shown  sche¬ 
matically  in  the  diagram  of  Fig.  4-9. 


Table  4-1 t 


Data  for  a  2nd  Order  Sub-harmonic  Oscillation  Response 
Curve,  in  which:  f(E)  *  2E  -  0.0048E3;  K  =  4,./-  0.25; 
and  9^(t )  «  2.5  Cos(  '  t  +  9) 


(a)  Frequency  Increasing: 


quency  of 

Amplitude  of 

Amplitude  o 

Input 

1st  Harmonic 

2nd  Sub-ham 

(  "  ) 

c  &,) 

( ^ ) 

0.628 

0.22 

/ 

0.754 

0.33 

/ 

0.880 

0.45 

/ 

1.000 

0.50 

/ 

1.130 

0.63 

/ 

1.260 

1.00 

/ 

1.380 

1.15 

/ 

1.510 

1.50 

/ 

1.640 

2.00 

/ 

1.760 

2.50 

/ 

1.890 

3.00 

/ 

1.950 

3.25 

/ 

1.980* 

10.50 

/ 

2.010 

10.00 

/ 

2.200 

9.00 

/ 

2.510 

7.30 

/ 

2.810 

6.00 

/ 

2.950 

/ 

/ 

3.140 

5,00 

/ 

3.460 

4.50 

/ 

3.770 

4.00 

/ 

4.090 

3.75 

/ 

4.40C 

3.50 

/ 

4.710 

3.25 

/ 

5.020 

/ 

4.75* 

5.150 

/ 

4.50 

5.270 

/ 

3.50 

5.340 

2.50 

/ 

Remarks 


*lst  harmonic 
upper  jump 


v2nd  sub-harmonic 
upper  jump 
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(b)  Frequency  Decreasing: 


Frequency  of 
Input 


Amplitude  of 
1st  Harmonic 


(  -  ) 

(  *  ) 

(  A 

5.340 

2.50 

/ 

5.270 

/ 

3.50 

3.150 

/ 

4.50 

5.020 

/ 

4.75 

4.710 

/ 

5.50 

4.400 

/ 

6.50 

4.  COO 

/ 

7.80 

3.770 

/ 

9.50 

3.460 

/ 

11.00 

3.140 

/ 

13.50 

2.950 

/ 

14.50* 

2.810 

6.00 

/ 

2.510 

7.25 

/ 

2.200 

8.75 

/ 

2.010 

10.00 

/ 

1.890 

11.00 

/ 

1.760 

12.25 

/ 

1 . 640 

13.00 

/ 

1.510 

15.00 

/ 

1.380 

17.00 

/ 

1.260 

21.20 

! 

1.130 

24.00 

/ 

1.000 

25.50 

/ 

0.880 

0,45* 

/ 

0.754 

0.33 

! 

0.628 

0-22 

/ 

Amplitude  of 

2nd  Sub-harmonic 
) 


Remarks 


*2nd  sub-harmoni< 
lower  jump 


*lst  Harmonic 
lower  jump 
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4*7  Analog  Computer  Analysis  for  a  Ird  Order  Stib-harrnonlc  Oscillation 
System: 

Consider  a  system  with  a  block  diagram  shown  in  Figo  4-10: 


r 


>  * 


Fig.  4-10  Block  Diagram  of  a  3rd  Sub-harmonic 
Oscillation  System 


in  which  assume  the  nonlinearity  of  system  is  an  "ON-OFF”  element,  a 
sketch  of  the  characteristic  for  the  computer  simulation  is  shown  in 
Fig.  4-11,  it  is  approximated  from  three  straight  lines;  one  of  which 
with  a  slope  of  15.0;  one  of  which  with  a  slope  of  0.1;  and  the  other 
is  zero: 

The  circuit  setting  up  in  the  computer  is  shown  in  Fig.  4-12,  all 
components  to  be  used  are  schematically  indicated  with  the  diagram.  The 
characteristic  curve  from  the  computer  is  shown  in  Fig.  4-13. 


r 


Fig.  4-12  Function  Generator  Diagram  for  an 
"ON-OFF"  Nonlinearity 
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I 


ig. 


4  ii  Cuurac tori stic  Sketch  of  Ah 
,!Cn-0i J'  Nonlinear  hlcncnt 


..  2 


I 
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An  operational  diagram  tor  setting  up  the  analog  computer  is  shown 
in  Fig.  4-14,  in  which  the  value  of  all  components  to  he  used  are 
schematically  indicated  in  the  diagram.  The  scaling  factor  for  the 
computer  setup  is  the  same  as  diagram  of  Fig.  3-20  to  be  used. 

The  results  from  the  computer  are  condensed  in  Table  4-2,  and  with 
a  response  curve  shown  in  Fig.  4-15.  It  is  shown  that,  there  is  a  pair 
of  jumps  for  each  harmonic  oscillation  response  curve,  these  jump  pheno- 
menas  within  frequency  range  is  accompanied  shown  with  the  diagram  of 
Fig.  4-15. 
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Table  4-2:  Data  for  a  3rd  Order  Sub-harmonic  Oscillation  Response 

Curve  from  an  MON-OFF"  Nonlinearity  System  and;  '  =  0.4 
K  =  10.0;  ©v^(t)  *  8.0  Cos(’ot  +  6) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Amp  lit due  of 

Amplitude  of 

Remarks 

Input 

1st  Harmonic 

2nd  Sub-har. 

3rd  Sub-har. 

(to  ) 

(  <=>  ) 

( 

( *$) 

0.628 

0.60 

/ 

/ 

0.880 

1.25 

/ 

/ 

1.130 

2.25 

/ 

/ 

1.380 

3.50 

/ 

/ 

1.640 

4.’5 

/ 

/ 

1.760 

5.75 

/ 

/ 

1.890 

7.50 

/ 

/ 

2.010* 

50.00 

/ 

/ 

*lst  Harmonic 

2.140 

46.00 

/ 

/ 

Upper  Jump 

’.270 

42.50 

/ 

/ 

2.390 

39.00 

/ 

/ 

2.510 

35.00 

/ 

/ 

2.630 

31.50 

/ 

/ 

2.760 

28.00 

/ 

/ 

2.890 

26.00 

/ 

/ 

3.010 

22.50 

/ 

/ 

3.140 

20.00 

/ 

/ 

3.260 

/* 

/ 

/ 

^Transient  and 

4.270 

/* 

/ 

/ 

unstable 

4.390° 

/ 

/ 

/ 

""Transient  from 

4.520 

/ 

/ 

/ 

1st  Harmonic 

4.650 

/ 

/ 

/ 

to  2nd  Sub-har. 

4.770 

/ 

/ 

/ 

4.990° 

/ 

/ 

/ 

5.020* 

/ 

/ 

/ 

5.140 

/ 

20.00 

/ 

*2nd  Sub-ha i. 

5.270 

/ 

18.00 

/ 

Upper  Jump 

5.400 

/ 

16.70 

/ 

5.520 

/ 

16.00 

/ 

5.650 

/ 

15.00 

/ 

5.770 

/ 

13.80 

/ 

5.900 

/ 

13.50 

/ 

6.030 

/ 

12.70 

/ 

6.150 

/ 

12.30 

/ 

6.280 

/ 

11.90 

/ 

6.400 

/ 

11.30 

/ 

6.530 

/ 

11.00 

/ 

6.650 

/ 

10.60 

/ 

6.780 

/ 

10.20 

/ 

6.910* 

/ 

10.00 

/ 

*3rd  Sub-i  r. 

7.030 

/ 

/ 

35.00 

Upper  Jump 

7.160 

/ 

/ 

32.50 

7.280 

/ 

/ 

28.50 

7.410 

/ 

/ 

27.20 

7.530 

/ 

/ 

26.00 
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(b)  Frequency  Decreasing: 


Frequency  of 
Input 
(u?  ) 

7.530 
7.410 

7.280 
7.160 
7.030 
6.910 
6.780 
6.650 

6.530 

6.400 

6.280 
6.150 
6.030 

5.900 

5.770 
5.650* 

5.520 

5.400 

5.270 

5.140 
5.020 

4.900 

4.770 
.  4.650 

4.520 

4.390 
4.27C  °* 

4.140 
3.260° 

3.140 
3.010 
2.890 

2.760 
2.630 

2.510 

2.390 

2.270 

2.140 
2.010 
1.830 

1.760 
1.640* 

1.510 
1.380 
1.130 
0.880 
0.628 


Amplitude  of 
1st  Harmonic 
(4  ) 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

20.00 

22.50 

26.00 

27.50 
31.00 
35.00 
39.00 
42.00 

45.50 
51.00 
55.00 

52.50 
5.00 

4.25 
3.50 
2.30 

1.25 
0.60 


Amplitude  of 
2nd  Sub-har. 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

13.75 

15.00 

16.00 

16.70 

17.50 

19.50 

21.50 
23.00 
25.00 
32.00 
36.00 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


Amplitude  of 
3rd  Sub-har. 

26.00 

27.20 

28.50 

30.50 

32.50 
35.00 

36.50 

37.50 
40.00 

42.50 
43.75 
47.00 
50.00 

52.50 
56.00 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


Pema  rks 


*3rd  Sub-har. 
Lower  Jump 


°*2nd  Sub-har. 
Lower  Jump 
and  Transient 
Unstable  Region 


*lst  Harmonic 
Lower  Jump 
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•f*C  6-0  7C 


CHAITEK  V 


FORCED  OSCILLATION  FOR  THE  HIGHER  ORDER 
FEEDBACK  CONTROL  SYSTEM 
5-1  General  Description; 

Up  to  the  Chapter  4*  only  forced  oscillations  m  the  case  of  the  2nd 
order  system  are  discussed.  In  the  extension  to  higher  order  systems* 
consider  a  block  diagram  shown  in  Fig.  5-1; 


o 


C;  N 


Cci* ) 


a. 


in  which  G(s) 
of  the  system* 


Fig.  5~1  Block  Diagram  of  a  Feedback 
Control  System 

is  an  equivalent  transfer  function  of  the  linear  element 
it  can  be  written  in  the  form: 

2M 


G(s) 


K 


F(5) 


(5*1) 


and  with  a  result  for  the  system: 

P(s)  E  +  Q(s)  Kf (£)  =  P(s)  (5-2) 

Where  P(s)  may  be  any  order  of  s;  as: 

P (s) ■  —  sn  +  b  sn  1  +  b0sn  2  + . +  b  (5*3) 

12  n 

Similarity  Q(s)  may  be 

,  ,  .  m  m-1  m-2  . _  . v 

Q(s)  s  s  +  c^s  4  C2®  + . o.+  cfn  (5*4) 

In  general  the  frequency  responsive  element  3(s)  following  the  non¬ 
linearity  is  of  such  a  nature  as  to  attenuate  the  high  frequency*  that 
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means  the  order  of  P(s)  should  bu  higher  than  the  order  of  Q(s).  The 
order  of  system  Is  n.  Inserting  equations  (5-3)  and  (5-4)  Into  (5-2)  and: 


.dl  _ 

i  e  7  1  .If 


aM ' 

lit  >  r 


rU  ) 


,»  d'-1 

Ot  f  *  di*-i  r  rb  1  )  (5-5) 

Equation  (5-5)  is  a  general  differential  equation  form  of  the 

forced  oscillation  of  any  order  feedback  control  system. 

5-2  Harmonic  Linearization  and  Iteration  Method  for  Investigating  the 
Equivalent  Gain  of  a  Nonlinearity: 

After  a  restoring  force  function  of  the  nonlinear  element  has  been 

defined,  the  equivalent  complex  gain,  N(E^)  of  the  nonlinear  element  is 

determined  by  the  ratio  of  the  complex  amplitude  of  the  first  harmonic 

of  the  output  to  the  amplitude  of  the  harmonically  varying  Input  quantity. 

Recall  the  block  diagram  of  Fig.  5-1,  if  the  system  is  a  2nd  order 

3 

system  and  with  a  restoring  forcing  function;  f(E)  ■  a^E  +  a^E  ,  and: 


K  N(E)  E  =  6 
s  (s  ) 

or: 

E  +  E  +  K  N(E)  E  =  9r  +  ^  &r 
Assume  the  input  function  is 

f~\  s  F  Cos(h  t  +  G) 
and  one  solution  of  E: 


E  =  Ej  Cos  uo  t 


and  with  a  result: 

(  ou-1  Ex  +  K  N(E1)E1)Costo  t  «  J  EjSinw  t 


-A  Cos  (a,  t  +  0) 


(5-6) 


(5-7) 


(5-8) 


(5-9) 


-  BSinOo  t  +  9)  (5-10) 

Comparing  the  result  with  the  equation  of  (3-42)  the  equivalent  gain: 


N(Ej) 


(5-11) 
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For  a  general  lor™  expression;  the  equivalent  gain  of  a  non¬ 
linearity  in  terms  of  the  amplitude  of  input  Is  devtlrped  by  Za 
Tsypkin.  The  result  for  a  symmetrical  nonlinearity  with  input  of  equation 
(5-9)  is: 

N(El)  *  [  f<V  +  f<Ei/2 >J  (5-12) 

If  more  accuracy  is  required,  the  equivalent  gain  N(Ep  may  be; 

N(E1}  "  3E^ - i  f(El)  4  f<El/2>  +  ,3  f(,3  Ej/2)  (5-13) 

For  the  case  of  a  nonlinearity  with  a  restoring  force  function* 

3 

f(E)  «  a^E  +  a^E  ,  the  results  from  equations  (5-12)  and  (5-13)  are  ex¬ 
actly  the  same  as  equation  (5-11).  Hence,  for  any  nonlinearity  the  case 
of  first  harmonic  input  the  equivalent  gain  can  be  expressed  either  in  the 
form  of  equation  (5-12)  or  (5-13). 

Equation  (5-12)  or  (5-13)  not  only  substantially  simplifies  the 
calculation  for  a  second  order  forced  oscillation,  it  is  also  probably 
a  most  important  for  investigating  self  oscillation  or  forced  oscillations 
in  a  higher  order  system  with  single  or  multiple  nonlinearities  in  a 
feedback  control  system. 

5-3  Forced  Oscillations  for  a  Higher  Order  Feedback  Control  System: 
Consider  a  system  with  a  linear  transfer  function: 

G(s)  a  _ K(s  4  V  (5-14) 

s(s  +  Pj)(s  +  P2)(s  +  P3) 


or  in  the  form: 


G(s) 


K(s  c) 


4  3  2 

s  -J-  b,  s  +  b^s  +  b^s 


(5-15) 


From  the  equation  (5-5),  the  differential  equation  for  the  system 


of  equation  (5-15): 

.4 


JL-.  .  K  _!.)£  +  K(  t  c)f(E) 

(  —  +  bl.  3  +  b2  2  4  b3  dr  dt 

dt 


dt 


dt' 
=  ( 


4  3  2 

d  dJ  l  ,  d__  ,  .  d  )  6^ 

dt4  1  dt3  2dt  3  dt 


(5-16) 


0  1 


in  which, 


E  •=  E jCos  i* /  L 
s  FCos(^t  +  G) 


and ; 


1  3 

f(E)  *  E^Cos  ot  +  £  (3Coscot  +  Cos3  *  t) 

Inserting  these  equations  into  equation  (5-16)  and  with  a  result! 


(&lCU+-y  fc'io  t  Kd'Li-',  r  of  r<Z,  C-  cf  )  ~ 

r  (  h.  L  ,  h  a  -  i,  £ri  CO  —  /\  4$  £*  0/.  bit a  v 

r  -j  Kdj  i-  3  t,c£>  rWfc  -  j-n'J  MJ  sic*  Wc 
~  hvV  [(A/  (' d  ^  ~h~)  u>£ft  r  bj  be 


-  fu.  /  it)  ( ci)  *— b  *■  X'tw  &  [b  b .  )  Cct  &  J 

(5-17) 

Neglect  the  higher  order  harmonic  terms  and  equate  the  coefficients 
of  Cos  a)  t  and  Sin  lit: 

£t,( CO'1  -  b  .  i*.  '  r  A  CL,  C  T  i  i  ^  1 

/  <•  c*'0  -t  (b.-  '■  (5-18) 

E,C bi  -{h  :  ,  t  jj-  Kd)  £,*)]& 

~  ?  •  •  '  !  w:  ~  lbtv)"-  f)c\^0)  (5-19) 


Squaring  equations  (5-18)  and  (5-19)  and  adding: 


£T<  -x 


it 


i  <~ 


VS1-  7i^  vv- 


—  j  v  /  iv'(i  r  ->  t  fb|li/  J 

(5-20) 

Equation  (5-20)  is  the  equation  of  response  curve  of  the  system  of 
equation  (5-15).  If  the  constants  of  system  are  know,  a  response  curve 
in  the  plane  can  be  plotted. 

For  the  investigation  of  conditions  for  existence  of  forced  oscil¬ 
lation  rearrange  the  equation  (5-20)  and  put: 
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X 


=  a. 


Hence : 


r-r  .  ft 


'V  ; 


(5-21) 


or; 


V  T 


‘ 


Z*/  Js,~ 

t  [  L*  I  !\  ±  -r  n'  /  *  r  -  A  s  (5-22) 

r 

The  condition  for  existence  of  forced  oscillation  is  found  from 
equation  (5-22)  by  means  that  the  solution  of  u  should  be  real.  It  is 
to  be  noted  that  the  value  of  X  is  the  equivalent  gain  of  nonlinearity, 
its  value  should  be  larger  than  zero,  as  it  was  proved  in  Chapter  3. 

For  the  phase  investigating  recall  the  equation  (5-19)  and  re¬ 
write  : 

v  (&  *  -h .  )  '  -  i 

1  i  •  ,  |  ■ 

— -  '  i.  '<  ~t-;  i\  4f  >  y  -i  '  (5-?3) 

or: 


<5,n  C  (V  -J  )  -  — - 

/  L  u 


j  •  .  <  c  ^  /- 

*  '■  - /  t-  *r  r  i,  .  •  • 


(5-2'*) 


hence : 


(9 


c 

~  /  ~ (  -V  ir ,  )  r  i  i,  , . 


r  y 


(5-2S) 


in  which  the  value  of; 

ip  —  Iuk  !Li  b.lX*'  k,)  .<  (  h,  1  .■/  (5-/-0) 

Equation  (5-25)  is  shown  that  the  value  of  8  is  a  function  of  «. 
and  the  amplitude  of  first  harmonic  of  error  signal,  if  the  system 
constants  is  fixed. 
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If  we  use  the  equivalent  gain  for  the  system  investigation,  and 


with  a  result : 


4 


U  1  y  v  t,  ) 


(5-27) 


/  t  ME,  )  iL  ) 


or : 


a 

& 


4  c  )  \ 


(5-28) 


Where  is  the  amplitude  of  first  harmonic  of  error  signal  of  forced 
oscillation. 

Equation  (5-28)  is  similar  to  a  linear  system,  the  techniques  for 
linear  system  are  applicable. 

5-4.  Analog  Computer  Analysis  for  a  Higher  Order  System: 

Consider  a  feedback  control  system  with  a  transfer  function  and  a 
nonlinearity  are  shown  in  Fig.  5-2: 


Fig.  5-2  Block  Diagram  for  a  Typical 
Feedback  Control  System 


in  which  the  restoring  force  function  is: 


f (E)  =  2E  +  0.020E3 


(5-29) 


The  characteristic  curve  of  restoring  force  function  of  nonlinear¬ 
ity  from  the  calculation  is  shown  in  Fig.  5-3.  Where: 


E  =  1.00;  2.00;  3.00;  4.00;  5.00;  6.00;  7.00;  7.50; 


f (E)  =  1.98;  3.84;  5.46;  6.72;  7.50;  7.68;  7.15;  6.60: 
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I 


[ 


Fig. 5- 3  Characteristic  Curve  of 
f(E)  ~2B  ~  0 . 0201i-?  Fret* 
Calculation 


$5 


For  the  simulation  ot  the  restoring  force  function  of  nonlinearity 
from  the  Computer,  it  is  approximated  by  three  straight  lines,  one  of 
which  with  a  slope  of  2;  one  of  which  with  a  slope  of  1;  and  the  other 
with  a  slope  of  zero.  The  circuit  setup  in  the  computer  is  shown  In 
Fig.  5-4,  all  components  used  are  schematically  in  the  diagram. 

The  characteristic  curve  of  restoring  force  function  from  the  com¬ 
puter  is  shown  in  Fig,  5-5;  (a)  shows  the  waveforms  of  the  input  and  the 
output;  (b)  shows  the  characteristic  curve  of  the  nonlinearity. 


U  .ti'  )  /  " 

~  ^  v*  -  n 

n  <  i 

t 

i 

n,  )\  X 

.  •  E;, 

V-  ^-j 

/T,  yh<) 
wVva**- 


Fig.  5~4  Analog  Computer  Setup  for 

f (E)  «  2E  -  0.020E3  Simulation 

The  operational  block  diagram  and  the  operational  circuit  diagram 
are  shown  in  Fig.  5-6  and  Fig.  5-7.  The  all  components  used  in  the 
operational  circuit  diagram  are  seehemat Ically  in  the  diagram. 

The  results  from  the  computer  are  condensed  in  Table  5-1,  and  with 
response  curves  shown  in  the  plane  and  plane  are  show  in  Fig.  5-S 
and  Fig.  5-9  respectively. 
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Fig.  5-6  Operational  Block  Diagram  setup 
in  the  Analog  Computer 


98 


Table  5-1.  Data  for  Response  Cut^e  from  Computer  for  a  Highet  Order 

System;  in  which  f (E)  -  2E  -  0.020E3,  6/  (t)  -  3.0  Cos<Vt  +6) 

(a)  Frequency  Increasing: 


Frequency  of 

Amplitude  of 

Phase  of 

Amplitude 

of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

) 

(  £?,) 

(  4'  ) 

(  <p  ) 

(  C  ) 

0.314 

0.40 

270.00 

3.10 

0.00 

0.377 

0.55 

270.00 

3.30 

0.00 

0.440 

0.75 

262.00 

3.40 

0.00 

0.503 

0.90 

255.00 

3.50 

5.80 

0.565 

1.10 

246.00 

3.60 

9.00 

0.628 

1.25 

241.00 

3.80 

12.50 

0.755 

1.70 

235.00 

4.30 

20.00 

0.880 

2.50 

220.00 

5.00 

24.50 

1.000 

7.40 

105.00 

6.30 

107.50 

1.130 

6.90 

64.00 

5.10 

120.00 

1.260 

6.50 

37.50 

4.25 

137.50 

1.380 

6.10 

25.00 

3.50 

155.00 

1.510 

5.60 

15.00 

3.00 

165.00 

1.640 

5.25 

10.00 

2.50 

172.00 

1.760 

4.90 

7.50 

2.00 

180.00 

1.890 

4.50 

5.00 

1.70 

186.00 

2.200 

4.00 

0.00 

1.20 

195.00 

2.500 

3.50 

0.00 

0.85 

210.00 

2.810 

3.30 

0.00 

0.60 

220.00 

3.140 

3.20 

0.00 

0.40 

232.00 

3.460 

3. 10 

0.00 

0.30 

245.00 

3.770 

3. 00 

0.00 

0.25 

260.00 

4.090 

3.00 

0.00 

0.20 

270.00 

(b)  Frequency 

Decreasing 

Frequency  jf 

Amplitude  of 

Phase  of 

Amplitude 

of 

Phase  of 

Input 

Error  Signal 

Error 

Output 

Output 

(  »  ) 

(  *,) 

(  ^  ) 

(  ) 

(  6  ) 

2.510 

3.50 

0.00 

0.78 

210.00 

2.200 

4.00 

0.00 

1.20 

195.00 

1.890 

4.50 

5.00 

1.65 

185.00 

1.760 

4.90 

7.50 

2.00 

180.00 

1.640 

5.25 

10.00 

2.50 

171.00 

1.510 

5.70 

15.00 

3.00 

162.00 

1.380 

6.10 

25.00 

3.50 

155.00 

1.260 

6.60 

37.50 

4.25 

137.00 

1.130 

6.60 

64.00 

5.00 

121.00 

1.000 

7.40 

105.00 

6.25 

106.00 

0.943 

/ 

/ 

6.80 

78.00 

0.880 

6.70 

150.00 

6.00 

25.00 

0.754 

1.80 

235.00 

4.30 

20.00 

0.628 

1.25 

248.00 

3.80 

16.00 
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CHAPTEF  VI 


CONCLUSIONS 

The  atm  of  this  dissertation  was  to  present  a  method  for  investigat¬ 
ing  the  conditions  for  existence  of  the  forced  oscillations  of  a  nonlinear 
feedback  control  system.  This  method  is  called  ’’harmonic  linearization 
method”,  it  is  based  on  the  harmonic  balance  and  iteration  method.  From 
this  method,  the  conditions  for  existence  of  the  fundamental  frequency 
and  the  sub-harmonic  of  order  2  and  3  are  investigated,  and  also  a  equi¬ 
valent  gain  of  the  nonlinearity  in  terms  of  input  amplitude  and  frequency 
has  been  developed. 

The  term  of  the  ’’response  curve”  as  used  here  is  defined  as  the  combina¬ 
tion  of  two  components,  one  is  the  attenuation  frequency;  it  is  either 
shown  in  the  ’’Error  vs  Frequency”  plane  or  ’’Output  vs  Frequency”  plane, 
and  the  other  is  the  phase  frequency  response,  it  is  with  a  difference 
value  from  the  different  representation  of  attenuation  frequency  response. 

The  results  of  this  analysis  show  that  a  nonlinear  feedback  control 
system  may  possess  more  than  one  type  of  forced  oscillations,  which  type 
will  exist,  it  is  depending  on  the  characteristic  of  the  system  and  the 
characteristic  of  input  forcing  function. 

For  any  nonlinear  feedback  control  system,  the  response  curve  of 
either  fundamental  or  sub-harmonic  forced  oscillations,  it  is  usually  a 
sudden  jump  in  amplitude  and  phase  when  a  jump  transient  take  place.  The 
range  of  jump  frequency  depends  on  the  value  of  damping  of  the  system  and 
the  amplitude  of  the  input;  usually  it  is  increased  as  the  damping  decreas¬ 
ing,  and  decreased  as  the  input  amplitude  decreasing.  A  vertical  tangents 
locus  investigating  for  calculating  the  jump  frequencies  in  terms  of  system 
constants  and  the  amplitude  of  input  was  investigated. 
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Due  to  the  resonance  phenomena  and  the  linear  approximations  of  re¬ 
storing  force  function,  there  is  some  difference  between  the  theoretical 
and  the  experimental  response  curves.  It  if  only  existant  in  the  jump 
frequency  range,  other  than  this  range  of  frequencies,  it  is  very  close 
to  each  other. 

Sub-harmonic  oscillations  are  undesired  for  a  nonlinear  feedback 
control  system,  it  is  actually  an  unstable  phenomenon.  Under  suitable 
conditions,  a  multiple  order  sub-harmonic  forced  oscillation  may  exist 
in  a  system  with  a  constant  input  forcing  function.  There  Is  also  a 
jump  phenomena  for  each  order  of  sub-harmonic  oscillations.  The  condi¬ 
tions  for  the  existence  of  the  2nd  and  3rd  order  sub-harmonic  forced 
oscillations  are  investigated,  and  also  theoretical  response  curve  equa¬ 
tions  for  each  case  are  developed. 

It  is  possible  an  important  advantage  that  the  design  and  analysis 
of  a  nonlinear  feedback  control  system  by  using  the  harmonic  lineariza¬ 
tion  method,  by  means  for  investigating  an  equivalent  gain  of  the  non¬ 
linearity.  In  this  case,  all  principles  and  techniques  applied  for  linear 
system  are  applicable. 
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